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^This  roport  doacrlbot  the  results  of  efforts  to  evaluate  sad  develop 
BMthods  of  flev  detection  In  laminated  fiber  composite  aerovehlde  structures 
using  Infrared  radiation  detactlon  techniques. 

While  not  exhibiting  the  quantitative  detection  capabilities  of  other 
nondestructive  evaluation  (NDE)  methods  such  as  ultrasonics  or  neutron 
radiography,  infrared  thermography  promises  to  provide  a  rapid  method  of 
damage  detection  in  large  structural  surfaces.  Two  methods  are  described:  \ 
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extenielly  applied  thermal  field  (EATF)  flaw  detection  where  surface  CHBot 
spots”’ due  to  subsurface  flaws  ara  craatad  by  baaciag  with  a  radiant  hast 
•ourca;  and  stress-generated  tbarnal  fiald  (SGTF)  flaw  dataetion  whara 
"hot  spots"  in  vlscoaiastie  composites  ara  craatad  at  straws  concentrations 
naar  structural  flaws  during  high  frequency,  low  nagnitude  cyclic  loading. 

Analytical  and  experimental  studies  ware  performed  on  8-  to  32-ply 
graphite/epoxy  leal nates  and  8-  to  16-ply  S2-glass/epoxy  laminates  containing 
delaminations  and  partial  through-holes  to  evaluate  and  develop  EATF  thermo¬ 
graphy.  Insults  show  that  the  conduction  method  of  heat  application  is 
unsatisfactory  for  NDB  purposes,  and  that  a  radiant  energy  source  in  excess 
of  one  solar  constant  will  be  required  for  bast  results. ^Flawe  ere  detect-  " 
able  by\this  method  up  to  7/8  through  the  thickness  of  32  ply  graphlte/epoxy 
and  16  ply  $2-glass/epoxy  laminates,  and  it  does  not  appear  that  limits  of 
flaw  detectlba  have  yet  been  found.  Surface  finish  (paint  thickness, 
abrasions,  ana'soll)  have  been  found  to  affect  results  in  certain  cases,  and 
methods  of  ellmlnatlsg  their  undesirable  effects  eta  presented. 

SGTF  evaluation  was  conducted  On  both  graphite/ epoxy  and  S2-glasa/epoxy 
laminates  using  a  closely  coordinated  program  of  analysis  and  experiment.  A 
semi-empirical  method  of  determining  heat  generated  in  eompoaltes  under 
cyclic  loading  has  been  developed,  and  reeults  indicate  that  the  nagnitude 
of  heat  generated  by  viscoelastic  materiel  behavior  may  be  marginal  for 
SGTF  flaw  detection.  Beat  generated  by  frictional  rubbing  of  crack  or 
delamination  surfaces  has  been  shown  to  be  several  orders  of  nagnitude 
larger  than  heat  due  to  material  viscoelasticity,  and  ultrasonic  excitation 
methods  suggested  by  other  researchers  My  prove  to  be  the  most  fruitful 
SGTF  technique.  j 


UNCLASSIFIED 

tacuwrY  ct,*swwc*Tww  os  thu  eessonw  i 


PREFACE 


The  research  reported  in  this  document  was  conducted  under  contract 
Mo.  N68335-79-C-1084.  This  research  is  a  continuation  of  that  begun  under 
contract  Mo.  N68335-78-M-5337,  AIRTASK  A34OOOOO/051B/8F414614O6,  to 
develop  thermographic  methods  of  rapid  flaw  detection  in  composite  aero- 
structures.  Dr.  P.V.  McLaughlin,  Jr.,  was  the  principal  Investigator,  Dr. 
E.V.  McAssey,  Jr,  was  co-investigator,  and  Messrs.  V.R.  Emany,  D.N.  Koert, 
and  M.  Spitzer  were  research  assistants.  All  authors  are  affiliated  with 
Vlllanova  University,  where  the  research  was  conducted.  Messrs.  R.C. 
Deltrich  and  P.V.  Clekurs  were  Naval  Air  Engineering  Center  technical 
monitors. 

The  authors  wish  to  thank  Dr.  John  F.  Kay  of  Owens -Corning  Fiberglass 
Corporation  for  providing  specially  fabricated  S2-glass/epoxy  test  panels 
for  use  in  this  program.  Thanks  are  also  due  to  Mr.  Alan  J.  Hannibal  of 
Lord  Corporation  for  providing  samples  of  graphite  and  Kevlar  composites. 
The  donation  of  test  machine  parts  to  Vlllanova  University  by  Mr.  Fred 
Giordano,  President,  KGM  Machine  and  Tool  Co.,  Inc.,  is  gratefully 
acknowledged. 


NAEC-92-157 


This  page  left  blank 
Intentionally. 


i 


/ 

i 


NAEC-92-157 


TABLE  OF  CONTENTS 

Section  Title _  Page 

PREFACE  .  1 

LIST  OF  ILLUSTRATIONS  .  4 

LIST  OF  TABLES .  7 

I  INTRODUCTION  .  9 

II  EXTERNALLY  APPLIED  THERMAL  FIELD  METHOD  ...  11 

A.  OBJECTIVE  AND  OUTLINE  OF  THE  EATF  STUDY  11 

B.  HEAT  TRANSFER  ANALYSIS .  11 

C.  EVALUATION  OF  SURFACE  FINISH  EFFECTS  .  .  12 

D.  THERMAL  FIELD  APPLICATION  METHODS  ....  17 

E.  EVALUATION  OF  FLAW  DETECTION  CAPABILITIES  22 

III  STRESS-GENERATED  THERMAL  FIELD  METHOD  ....  28 

A.  OBJECTIVES  AND  OUTLINE  OF  SGTF  STUDY  .  .  28 

B.  ANALYSIS  OF  HEAT  GENERATION  IN  CYCLICALLY 

LOADED  COMPOSITES .  28 

C.  EVALUATION  OF  CYCLICALLY  LOADED  GRAPHITE/ 

AND  GLASS/EPOXY .  43 

D.  FEASIBILITY  ASSESSMENT  OF  THE  SGTF  METHOD  47 

E.  ULTRASONIC  STRESS  GENERATION  METHODS  .  .  48 

IV  SUMMARY  AND  CONCLUSIONS .  50 

V  CURRENT  STATUS  AND  RECOMMENDATIONS .  55 

VI  REFERENCES .  56 


3 


NAEC-92-157 


LIST  OF  ILLUSTRATIONS 

Figure 

No.  _ Title _  Page 

1  THERMOVIDEOGRAPH  OF  BLACK/WHITE  NAVY  AVIATION  PAINT 

ON  ALUMINUM  PLATE .  65 

2  THERMOVIDEOGRAPHS  OF  WHITE/  ORANGE  PAINT  ON  ALUMINUM  .  66 

3  THERMOVIDEOGRAPHS  OF  WHITE  PAINT  ON  ALUMINUM  ....  67 

4  THERMOVIDEOGRAPHS  OF  COMPOSITE  PLATES  WITH  SINGLE  AND 

DOUBLE  COATS  OF  WHITE  NAVY  AVIATION  PAINT  .  68 

5  THERMOVIDEOGRAPHS  OF  SINGLE  COATS  OF  BLACK  AND  WHITE 

NAVY  AVIATION  PAINT  ON  COMPOSITE  PLATES  .  69 

6  THERMO VIDEOGRAPHS  OF  SINGLE  COATS  OF  ORANGE  AND  WHITE 

NAVY  AVIATION  PAINT  ON  COMPOSITE  PLATES  .  70 

7  THERMOVIDEOGRAPHS  OF  SINGLE  COAT  OF  WHITE  NAVY 

AVIATION  PAINT  ON  COMPOSITE  PLATES  .  71 

8  THERMOVIDEOGRAPHS  OF  SAMPLE  GR-D  WITH  CLEAN  ABRASIONS 

ON  WHITE  PAINT .  72 

9  SAMPLE  GR-D  WITH  CLEAN  ABRASION  AND  ABRASION  COVERED 

WITH  COAT  OF  CLEAN  MOTOR  OIL .  73 

10  GR-D  WHITE  PAINTED  SAMPLE  CONTAINING  CLEAN  ABRASION, 

ABRASION  WITH  HEAVY  COAX  OF  SOILED  GREASE,  AND  LIGHT 

COAT  OF  SOILED  GREASE  ON  GLOSSY  SURFACE .  74 

11  THERMOVIDEOGRAPH  OF  SAMPLE  CONTAINING  ABRASIONS  AND 

SOILED  GREASE  ON  NAVY  AVIATION  PAINT .  75 

12  SCHEMATIC  OF  UNIDIRECTIONAL  GRAPHITE/ EPOXY  PLATE 
CONTAINING  0.65-ca-DIA  PARTIAL-THROUGH-HOLE  USED 

IN  LOCAL  FIELD  APPLICATION  METHOD  EVALUATION  ....  76 

13  GRAPHITE/EPOXY  SAMPLE  OF  FIGURE  12  HEATED  WITH  HEAT 

GUN .  77 

14  GRAPHITE/EPOXY  SAMPLE  OF  FIGURE  12  WITH  HEAT  APPLIED 

BY  SURFACE  STRIP  HEATERS  .  78 

15  THEORETICAL  TEMPERATURE  PROFILE  FOR  GRAPHITE/ EPOXY 

SAMPLE .  79 

16  PHOTOGRAPH  OF  1320  W  SPACE  HEATER .  80 

17  EATF  SAMPLE  MOUNTING  AND  TEST  FIXTURES .  81 

18  THEORETICAL  SURFACE  TEMPERATURE  RESULTS  OF  25.4  AND 

12.7  mm  SQUARE  DELAMINATIONS  AT  MID-PLANE  OF  32-PLY 
GRAPHITE/EPOXY  PANEL  .  82 

19  RADIANT  HEAT  SOURCE  DESIGN  GEOMETRY  .  83 


4 


NAEC-92-157 


LIST  OF  ILLUSTRATIONS  (CONT'D) 


Figure 

No. 

Title 

Page 

20 

SOURCE  TEMPERATURE  versus  WIRE  DIAMETER  FOR  SINGLE 

SOURCE  RADIANT  HEATER  . 

84  ! 

21 

FINAL  DESIGN  GEOMETRY  FOR  HIGH  POWERED  HEATER  . 

85 

22 

RADIANT  HEAT  SOURCE  DESIGNED  AND  BUILT  FOR  EATF 

STUDIES  . 

86  1 

23 

VERTICAL  MOUNTING  ARRANGEMENTS  OF  GRAPHITE/ EPOXY 

PANEL  WITH  DEFECTS . 

87  \\ 

24 

INFRARED  CAMERA  ARRANGEMENT  AND  GR/EP  TEST  PANEL  FLAW 
DESIGNATION  . 

r 

88  !' 

25 

THERMOVIDEOGRAPH  OF  8-PLY  GRAPHITE/ EPOXY  PANEL  CON¬ 
TAINING  DELAMINATIONS  1  PLY  AND  4  PLIES  DOWN  FROM 

THE  SURFACE . . 

i 

89 

26 

THERMOVIDEOGRAPH  OF  FLAWS  B2  AND  C2  IN  32-PLY  GRAPHITE/ 
EPOXY  PANEL,  60  SECONDS  AFTER  COOLING  . 

90 

27 

THERMOVIDEOGRAPHS  OF  FLAWS  A2  AND  B2  IN  32-PLY  GRAPHITE/ 
EPOXY  PANEL  . 

1 

28 

THERMOVIDEOGRAPHS  OF  FLAWS  B1  AND  Cl  LOCATED  4  PLIES 

DOWN  IN  32-PLY  GRAPHITE/SPOXY  SAMPLE . 

1  • 

92  ; 

29 

THERMOVIDEOGRAPHS  OF  FLAWS  A2,  B2,  AND  C2  LOCATED  HALF 
THROUGH  32-PLY  GRAPHITE/EPOXY  PANEL  . 

i 

93 

30 

THERMOVIDEOGRAPHS  OF  FLAWS  A1  AND  Cl  LOCATED  28  PLIES 
DOWN  IN  32-PLY  GRAPHITE/ EPOXY  PANEL  . 

■: 

94 

31 

THERMOVIDEOGRAPHS  OF  FLAWS  B1  AND  Cl  LOCATED  28  PLIES 
DOWN  FROM  FRONT  SURFACE  IN  32-PLY  GRAPHITE/ EPOXY  PANEL. 

95 

32 

DEFECT  MAP  OF  ANISOTROPIC  32-PLY  [02/.-H5]8  GR/EP  PANEL 

AS  SEEN  FROM  BLIND  SIDE  . 

96 

33 

THERMOVIDEOGRAPHS  OF  BLIND  SIDE  OF  ANISOTROPIC  [02/±45]s 
GR/EP  PANEL  WITH  PARTIAL-THROUGH-HOLES  AND  DELAMINA¬ 
TIONS  . 

97 

34 

DEFECT  MAP  OF  16-PLY  [0/±45/90]s  S2-GLASS/EP0XY  PANEL 

AS  SEEN  FRGM  BLIND  SIDE  . 

98 

35 

THERMOVIDEOGRAPHS  OF  PARTIAL-THROUGH-HOLES  (PTH)  AND 
DELAMINATIONS  (DELAM)  IN  16-PLY  [0/±45/90]8  S2  GLASS/ 
EPOXY  PANEL  . 

99 

36 

PROCEDURE  FOR  CALCULATING  THE  TEMPERATURE  DISTRIBUTION 

IN  A  COMPOSITE  LAMINATED  STRUCTURE  UNDER  CYCLIC 

LOADING  . 

100 

5 


NABC-92-157 


LIST  OF  ILLUSTRATIONS  (CONT’D) 

Figure 

No.  _ _ Title _  Page 

37  DETERMINISTIC  PHASE  GEOMETRIES  FOR  UNIDIRECTIONAL  FIBER  101 

COMPOSITES . . 

38  UNIDIRECTIONAL  FIBER  COMPOSITE  COORDINATE  SYSTEM  AND 

AVERAGE  LAYER  STRESSES .  102 

39  METHOD  OF  DETERMINING  PHASE  STRESSES  DUE  TO  LAYER 
TRANSVERSE  NORMAL  STRESS  USING  SUPERPOSITION  OF 
TRANSVERSE  SHEAR,  PLANE  ISOTROPIC  STRAINING,  AND 

AXIAL  COMPRESSIVE  STRAINING  .  103 

40  STRESS-GENERATED  THERMAL  FIELD  TEST  APPARATUS .  104 

41  TYPICAL  HYSTERESIS  LOOPS  OBTAINED  DURING  ENERGY  DISSI¬ 

PATION  CALIBRATION  TESTS  OF  GRAPHITE/ EPOXY  AND  S2- 
GLASS/EPOXY  COMPOSITES  .  105 

42  ENERGY  DISSIPATED  PER  CYCLE  -  CELION  6000/NARMC0  5213 

GRAPHITE/EPOXY  COMPOSITE .  106 

43  ENERGY  DISSIPATED  PER  CYCLE  -  3M  SP-250-S2  S 2-GLASS/ 

EPOXY  COMPOSITE .  107 

44  FINITE  ELEMENT  GRID  FOR  25.4-mm-WIDE  SAMPLE  CONTAINING 

6.35-nm-DIA  THROUGH-HOLE .  108 

45  a.  FINITE  ELEMENT  RESULTS  -  HORIZONTAL  NORMAL  STRESS, 

ox .  109 

b.  FINITE  ELEMENT  RESULTS  -  VERTICAL  NORMAL  STRESS,  oy.  110 

c.  FINITE  ELEMENT  RESULTS  -  SHEAR  STRESS,  T^y .  Ill 

46  a.  HEATING  RATE  DENSITY  ANALYSIS  RESULTS  FOR  GRAPHITE/ 

EPOXY .  112 

b.  HEATING  RATE  DENSITY  ANALYSIS  RESULTS  FOR  S 2-GLASS/ 

EPOXY .  113 

47  NODAL  PATTERN  USED  IN  NUMERICAL  HEAT  TRANSFER  ANALYSIS 

FOR  CYCLED  GR/EP  AND  S2-GL/EP  SAMPLES  .  114 

48  THERMOVIDEOGRAPHS  OF  [0/145/90]  COMPOSITE  SPECIMENS  WITH 

6-tnm  THROUGH-HOLE  UNDER  CYCLIC  TENSION .  115 

49  THERMO VIDEOGRAPHS  OF  [0/±45/90]  COMPOSITE  SPECIMENS 
CONTAINING  THROUGH-HOLES  AND  DELAMINATIONS,  UNDER 

CYCLIC  TENSION  .  116 

50  THERMOVIDEOGRAPHS  OF  [0]  GRAPHITE/EPOXY  SAMPLE  UNDER 

CYCLIC  TENSION  SHOWING  FRICTIONAL  HEAT  .  117 


6 


NAEC-92-157 


LIST  OF  TABLES 

Table 

No.  _ Title _  Page 

1  AGA  THERMOVISION  SYSTEM  680/102B  COLOR  SCREEN 

TEMPERATURES  AND  SENSITIVITIES  .  57 

2  ALUMINUM  SURFACE  FINISH  SAMPLES .  58 

3  FIBER  COMPOSITE  SURFACE  FINISH  SAMPLES  .  59 

4  CONSTANTS  FOR  USE  IN  DETERMINING  MATRIX  STRESSES 
DUE  TO  TRANSVERSE  NORMAL  COMPOSITE  STRESS, 

EQ.  (6) .  60 

5  MODEL  CALIBRATION  SPECIMEN  DIMENSIONS  AND 

MATERIAL  PROPERTIES . 61 

6  SGTF  CALIBRATION  AND  CALCULATION  RESULTS  -[±45] 

GRAPHITE/EPOXY  AND  S 2-GLASS /EPOXY . f  62 

7  HEAT  ENERGY  DISSIPATION  PROPERTIES  OF  GRAPHITE/ 

EPOXY  AND  S2-GLASS/ EPOXY  COMPOSITES  .  63 

8  COMPOSITE  THERMOELASTIC  PROPERTIES  USED  IN 
CALCULATING  SURFACE  TEMPERATURES  OF  [0/±45/90] 

SPECIMENS  WITH  6-nnn  THROUGH-HOLES .  64 


This  page  left  blank 
intentionally . 


NAEC-92-157 


I.  INTRODUCTION 

A.  One  of  Che  Important  missions  of  serovehlcle  ground  support  personnel 
Is  to  locate  structural  damage  sites  or  flaws  which,  if  not  repaired,  would 
in  time  lead  to  structural  failure.  Several  nondestructive  evaluation  (NDE) 
methods  are  currently  available  or  are  under  development  for  this  purpose; 
among  them  are  X-ray,  neutron  radiography,  sonic  or  ultrasonic  transmission 
or  reflection,  eddy  current,  etc.  Each  of  these  methods  currently  requires 
sophisticated  equipment  and  highly  skilled  technicians  for  operation  and 
data  interpretation.  Several  methods  (for  example,  C-scan)  are  not  at  this 
time  suited  for  examining  large  areas  of  aircraft  components,  and  other 
methods  also  have  their  limitations. 

B.  The  use  of  fiber  composite  materials  in  present  and  future  generation 
Navy  aircraft  requires  that  NDE  techniques  be  developed  to  handle  aero- 
structures  made  of  composites  as  well  as  conventional  metals  such  as 
aluminum  and  titanium.  This  poses  specific  problems,  since  composite 
structures  fail  in  a  different  manner  than  metal  structures.  For  example, 
impact  damage  to  composite  skins  may  consist  of  delaminations  between 
plies,  or  subsurface  shattering  of  the  composite  laminate;  and  both 
generally  occur  without  visible  damage  to  the  top  surface. 

C.  A  relatively  new  nondestructive  evaluation  method  has  undergone 
extensive  examination  at  Vlllanova  University  laboratories  following  the 
initial  feasibility  study  presented  in  reference  (a)*.  This  method  consists 
of  using  an  infrared  camera  to  identify  perturbations  in  structural 
surface  temperatures  that  may  be  caused  by  the  presence  of  a  defect  or 
damage  zone.  The  surface  temperature  differences  can  be  caused  by 
applying  heat  externally  with  a  radiant  heat  source  or  hot  air  heat  gun 
(externally  applied  thermal  field  (EATF))  or  by  applying  multiple  cycles  of 
loading  to  a  part  and  allowing  viscoelastic  behavior  or  internal  friction 

to  cause  "hot  spots"  in  the  vicinity  of  stress  concentrations  caused  by 
flaws  (stress  generated  thermal  field  (SGTF)).  These  studies  have  shown 
that  NDE  by  thermal  detection  can  be  used  for  both  composite  and  metal 
structures,  and  has  the  potential  for  large  surface  area  examinations 
using  ground  support  personnel  who  are  not  highly  trained. 

*Please  see  the  list  of  references  presented  on  page  56  of  this  report. 
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D.  A  closely  coordinated  analytical/experimental  program  was  conducted 
to  evaluate  thermography  as  an  NOE  tool  for  Navy  aircraft.  For  the  ex¬ 
perimental  portion  of  the  program,  temperature  gradients  were  induced  in 
flawed  specimens  of  graphite/epoxy,  S2-glaas/ epoxy ,  Kevlar /epoxy  and 
aluminum  by  applying  heat  externally  (EATF)  or  by  stress-cycling  (SGTP). 
The  resulting  transient  surface  temperatures  were  monitored  by  an  AGA 
Thermovision  System  680/102B  infrared  camera  with  both  black  and  white 
and  color  Isotherm  video  readout.  The  color  television  screen,  which  was 
used  for  the  majority  of  the  studies,  breaks  a  given  range  of  temperature 
down  into  ten  subranges,  and  shows  each  1/10  of  the  total  temperature 
range  as  a  different  color.  For  example,  a  10°C  gradient  from  20°C  to 
30°C  would  show  up  on  the  color  screen  at  10°C  sensitivity  as  (see  Table 
1): 


TEMPERATURE 

TEMPERATURE 

(°C) 

COLOR 

(°C) 

COLOR 

Above  29 

White 

23-24 

Light  Blue 

28-29 

Yellow 

22-23 

Green 

27-28 

Orange 

21-22 

Blue  Green 

26-27 

Red 

20-21 

Dark  Blue 

25-26 

Lavender 

Belov  20 

Black 

24-25 

Purple 

As  a  result,  borders  between  colors  represent  isotherms  (for  example:  the 
border  between  orange  and  red  indicates  surface  temperatures  of  27°C) . 
Perturbations  in  otherwise  smooth  Isotherms,  or  "hot  spots"  (concentric 
closed  Isotherms),  Indicate  flaws  interrupting  otherwise  uniform  heat  flow 
or  uniform  temperature  fields.  Figures  1,  5,  9,  and  33  are  typical 

examples.  Exact  temperatures  are  not  necessary  for  identifying  defects  hy 
perturbed  Isotherms,  only  a  temperature  difference  from  the  surrounding 
region  need  be  identified. 

E.  This  report  first  describes  the  heat  transfer  analysis  which  was 
used  to  analytically  predict  specimen  surface  temperatures,  then  presents 
results  of  studies  to  evaluate  surface  finish  effects,  EATF  thermal  field 
application  methods,  and  EATF  flaw  detection  capabilities.  Next,  the  SGTF 
heat  generation  analysis  and  test  program  is  described,  and  a  feasibility 
assessment  of  SGTF  methods  (Including  ultrasonics  as  a  heat  generation  tool) 
is  made.  Conclusions  are  presented  for  both  EATF  and  SGTF  methods,  and 
recommendations  are  made  for  work  necessary  to  develop  thermographic  NOT 
into  a  workable  Navy  aerostructure  test  method. 

in 
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II.  EXTERNALLY  APPLIED  THERMAL  FIELD  METHOD 

A.  OBJECTIVE  AND  OUTLINE  OF  THE  EATF  STUDY.  The  externally  applied 
thermal  field  research  had  as  its  objectives  to  develop  and  evaluate  heat 
application  techniques,  to  generate  a  data  base  for  sizes  and  types  of 
flaws  that  can  be  detected  by  EATF  techniques,  to  make  a  determination  of 
the  effects  of  surface  finish,  and  to  present  preliminary  conclusions  con¬ 
cerning  equipment  procedures.  In  order  to  accomplish  these  objectives,  it 
was  necessary  to  develop  a  numerical  heat  transfer  analysis  procedure  and 
codify  it  on  the  digital  computer.  This  heat  transfer  analysis  is  briefly 
described  in  the  next  paragraph.  It  is  followed  by  the  surface  finish 
effects  results,  including  evaluation  of  thickness  and  color  of  Navy 
aviation  paint  on  aluminum  and  epoxy  matrix  composites,  and  the  effects  of 
abrasions  and  foreign  substances  on  thermovideograph  infrared  camera  out¬ 
put.  Next,  thermal  field  application  methods  are  discussed,  and  the  ad¬ 
vantages  and  disadvantages  of  a  heat  gun,  heat  laiqp,  local  strip  heaters, 
and  radiant  heat  sources  are  presented.  Finally,  an  evaluation  of  flaw 
detection  capabilities  for  graphite/epoxy  and  glass/epoxy  composite  lami¬ 
nates  is  made,  including  size,  type,  and  location  limits  to  the  EATF 
technique . 

B.  HEAT  TRANSFER  ANALYSIS.  Analytical  studies  were  required  to  support 
the  experimental  phases  of  the  present  study.  In  particular,  a  multi¬ 
dimensional  transient  computer  program  was  developed  to  calculate  the 
transient  thermal  response  of  samples  exposed  to  both  externally  applied 
and  stress  generated  thermal  fields.  The  transient  thermal  response  pro¬ 
gram  was  developed  using  a  finite  difference  model.  In  this  model,  the 
sample  is  divided  into  a  series  of  isothermal  nodes.  The  transient  response 
of  the  1^  node  is  given  by  the  following  expression: 
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where: 

-  temperature  of  ith  node  at  end  of  time  step  AS 

“  temperature  of  it*1  node  at  the  beginning  of  time  step 

-  heat  dissipated  in  the  1th  node 
m^  “  mass  of  ith  node 

■  specific  heat  of  ith  n0de 
kji  ■  thermal  coupling  between  and  ic^  node 

In  this  procedure,  the  convective  heat  transfer  to  the  ambient  environ¬ 
ment  is  included  as  a  thermal  coupling  to  a  known  tenperature .  The  above 
equation  is  solved  for  each  node  at  each  time  step  from  the  initial  time 
to  the  final  time.  The  program  adjusts  the  time  step  used  in  the  calcula¬ 
tion  to  insure  that  stability  is  maintained. 

C.  EVALUATION  OF  SURFACE  FINISH  EFFECTS 

1.  NAVAL  AVIATION  PAINT  ON  ALUMINUM 

a.  The  Naval  Air  Engineering  Center  (NAVAIRENGCEN)  provided 
aluminum  samples  with  Navy  aviation  paint  for  the  first  phase  of  the 
surface  finish  evaluation.  The  samples  are  10  cm.  x  30.5  cm.  aluminum 
substrates  which  were  painted,  one-half  one  color  and  one-half  another. 

Since  the  purpose  of  the  test  was  surface  finish,  no  defects  were  put  into 
the  aluminum.  Table  2  presents  the  various  paint  finishes  examined.  All 
the  finishes  were  glossy  polyurethane  paint  which  will  be  used  on  all  Navy 
aircraft  in  the  future.  Each  sample  was  suspended  vertically  and  heated 
with  a  hot-air  gun.  After  being  heated  uniformly,  the  sample  was  allowed 
to  cool  by  free  convection  and  infrared  thermovideograph  pictures  were 
taken.  To  check  the  effect  of  reflections  from  external  sources,  a  60-watt 
incandescent  bulb  was  positioned  1.2  metres  from  the  sample  and  at  an  angle 
of  approximately  15  degrees  from  the  camera  optical  axis. 

b.  The  surface  finish  tests  yielded  a  number  of  interesting  results 
Figure  1  presents  the  results  for  Sample  A  which  was  a  black/white  sample. 
Figure  la  shows  the  sample  in  the  hot  condition.  The  light  color  on  the 
left  is  the  black  finish,  which  is  hotter  than  the  white  portion.  The 
sensitivity  is  set  at  20°c  from  black  to  white.  The  unevenness  along 

the  Interface  between  the  color  is  attributed  to  paint  thickness  variation. 
Figure  lb  shows  the  same  sample  but  at  a  later  time  so  that  the  sample  is 
cooler  than  in  Figure  la.  The  sensitivity  for  this  figure  is  5°C. 


NAEC-92-157 


Figure  lc  shows  Sample  A  under  the  Influence  of  the  light  bulb,  whose 
image  appears  In  the  upper  part  of  the  figure. 

c.  Figure  2  presents  results  for  Sample  C  which  was  white/orange. 

In  Figure  2a,  the  effect  of  color  and  paint  thickness  can  be  seen.  The 
nonuniform  pattern  can  be  attributed  to  different  emissivities  or  uneven 
paint  thickness.  Figure  3a  shows  a  similar  pattern  for  a  white/white 
sample  (Sample  E).  Since  the  white/white  sample  has  uniform  emlssivity, 
it  would  appear  that  the  unevenness  in  Figures  2a  and  3a  is  due  to  paint 
thickness  variations.  This  is  reinforced  by  the  fact  that  the  horizontal 
streaking  patterns  are  similar.  In  Figure  lb,  the  effect  of  emlssivity 

is  more  obvious  because  the  sample  is  of  uniform  intensity  but  slightly 
different  shading. 

d.  The  effect  of  emlssivity  differences  between  the  different 
color  finishes  would  be  to  raise  or  lower  the  local  temperature  and, 
therefore,  change  the  shade  of  the  picture.  In  general,  emlssivity  is  in¬ 
dependent  of  color  and  depends  upon  the  chemical  composition  of  the  paint. 
However,  the  shading  differences  in  Figure  lb  are  equivalent  to  only  1°C 
to  2°C,  which  could  be  related  to  slight  color  effects.  Differences  in 
paint  thickness,  on  the  other  hand,  would  cause  local  hot  spots  due  to  the 
thermal  resistance  of  the  paint.  These  hot  spots  would  tend  to  follow 
thickness  variations  rather  than  color  changes. 

e.  The  effect  of  the  light  bulb  can  be  seen  in  Figures  lc,  2b,  and 
3b.  In  all  cases,  the  effect  is  very  pronounced  and  would  overshadow  all 
other  inputs.  The  sample  surface  which  is  glossy  tends  to  reflect  the  heat 
of  external  sources.  Since  the  finish  used  by  Navy  aircraft  will  be  glossy, 
the  technician  using  an  infrared  camera  for  NDT  will  be  required  to  reduce 
the  background  of  high  energy  external  sources.  Lower  energy  or  more 
diffuse  sources  such  as  a  hot  runway  are  not  expected  to  be  a  problem. 

2.  NAVY  AVIATION  PAINT  ON  EPOXY  MATRIX  COMPOSITES 


a.  The  appearance  of  streaks  in  infrared  camera  pictures  due  to 
uneven  paint  thickness  on  aluminum  samples  led  to  the  concern  that  paint 
thickness  variations  might  be  picked  up  by  Navy  technicians  as  structural 
flaws  in  composite  aircraft  during  inspection.  The  polymeric  Navy  avia¬ 
tion  paint  has  a  relatively  lov  conductivity,  as  do  all  plastics,  while 
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the  aluminum  sample  has  a  very  high  conductivity.  It  was  reasoned  that  the 
streaks  or  hot  spots  on  painted  aluminum  were  caused  by  heat  being  conducted 
away  from  the  surface  of  thinly  painted  regions  by  the  aluminum,  while  the 
thicker  paint  regions  remained  hotter  due  to  the  larger  thickness  of  low 
conductivity  material.  If  this  were  true,  then  the  spotting  or  streaking 
seen  on  painted  aluminum  would  not  be  seen  on  painted  epoxy  matrix  compo¬ 
sites,  since  the  conductivity  of  the  paint  and  the  transverse  (through-the- 
thickness)  conductivities  of  the  composites  are  nearly  equal. 

b.  To  test  this  hypothesis  and  to  see  if  paint  thickness  variations 
on  composite  materials  could  result  in  hot  spots  similar  to  those  generated 
by  structural  flaws,  12  plate  samples  (6  each  of  unidirectional  graphite/ 
epoxy  and  Kevlar  49/epoxy)  of  epoxy  matrix  composites,  each  approximately 

6  cm  wide  x  28  cm  long  x  0.2  cm  thick,  were  sent  to  NAVAIRENGCEN  facilities 
where  they  were  unevenly  painted  with  black,  white,  and  orange  Navy  aviation 
paint.  Single  thicknesses  and  double  thicknesses  of  the  paints  were  applied 
side-by-side  and  on  top  of  one  another  according  to  the  scheme  shown  in 
Table  3.  Each  sample  was  photographed  at  room  temperature  and  at  an 
elevated  temperature  after  being  heated  by  a  radiant  heat  source.  Like  the 
painted  aluminum  samples,  the  painted  composite  samples,  as  expected, 
exhibited  an  infrared  response  which  varied  with  color  due  to  paint  pigment 
emls8ivlty  differences.  However,  also  as  expected,  the  large  paint  thick¬ 
ness  variations  did  not  appear  to  be  detected  by  the  infrared  camera. 

c.  Figure  4  shows  infrared  camera  thermovldeographs  of  graphite/ 
epoxy  (sample  Gr-A)  and  Kevlar/epoxy  (sample  K-A)  plates  with  single  and 
double  thicknesses  of  white  paint.  Both  the  graphite  and  Kevlar  samples 
show  vertically  uniform  temperature  distribution  in  the  region  of  the  change 
from  single  to  double  paint  thickness.  Figure  5  shows  thermovldeographs 
for  single  adjacent  (nonuniform)  coats  of  white  and  black  paints.  There 

is  no  dl8cemable  difference  in  heat  emission  of  the  two  paints  at  room 
temperature  (Figure  Sa).  At  elevated  temperatures  created  by  heating  with 
a  radiant  heat  source,  however,  there  is  considerable  difference  in 
Thermovision  output  between  the  black  and  white  paints,  with  the  black 
paint  appearing  as  much  as  2°C  hotter  than  the  white  paint  (Figure 
5b-c) .  There  is  not  any  evidence,  however,  of  spotting  or  streaking  due 
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to  uneven  paint  thicknesses .  (Spots  in  Figure  Sb  are  caused  by  uneven 
heat  distribution- the  spotting  is  not  reproducible.)*  Figure  6  shows 
similar  results  for  orange/white  paint  combinations.  Note  that  the 
emissivity  effects  due  to  color  are  not  as  evident  between  the  orange 
and  white  paints  (a  few  tenths  of  a  degree  C)  as  they  are  between  black 
and  white  paints.  Figure  7  presents  thermovideographs  of  one  (uneven) 
coat  of  white  paint  on  Kevlar/and  graphite/epoxy  samples.  Note  that 
streaks  or  spots  do  not  exist  due  to  uneven  paint  thicknesses. 

3.  ABRADED  PAINTED  SURFACES 

a.  During  a  preliminary  study  (reference  (a)),  it  was  observed 
that  dull  or  non-reflecting  areas  in  an  otherwise  shiny  surface  could  be 
Identified  in  the  infrared  camera  video  picture.  Navy  aircraft  are  likely 
to  receive  non-damaging  surface  abrasions  from  support  or  flight  personnel 
or  as  the  result  of  a  mission,  etc.  It  was  therefore  deemed  important  to 
determine  the  effect  of  abrasions  on  infrared  camera  video  readout.  The 
graphite/ epoxy  with  a  single  coat  of  white  paint  (aample  Gr-D)  was  first 
cleaned,  then  abraded  on  two  2.5  cm  square  areas  with  No.  120  grit  sand¬ 
paper.  The  abrasions  were  complete;  no  gloss  was  left  inside  the  designa¬ 
ted  areas,  but  the  paint  was  not  penetrated  through  to  the  graphite/epoxy 
substrate. 

b.  Thermovideographs  of  the  abraded  surface  are  shown  in  Figure 
8.  At  room  temperature  (Figure  8a),  the  abrasions  appear  slightly  cooler 
(about  0.2°C)  than  the  surrounding  glossy  surface.  When  heated  one  minute 
by  a  radiant  heat  source  and  allowed  to  cool  for  30  seconds,  there  is  no 
visible  difference  between  abraded  and  glossy  painted  surfaces  at  5°C 
Infrared  sensitivity  (Figure  8b).  After  cooling  for  7  minutes,  however, 
the  same  sample  shows  the  abrasions  appearing  about  0.2°C  lower  in  temper¬ 
ature  than  the  surrounding  glossy  area  (Figure  8c).  Zt  is  noted  that 
0.2°C  is  very  small,  and  is  near  the  Instrument's  lower  limit  on  resolution. 

4.  FOREIGN  SUBSTANCES  ON  FAINTED  SURFACES 

a.  In  addition  to  abrasions,  it  is  likely  that  portions  of  an  air¬ 
craft  surface  may  contain  foreign  substances  such  as  lubricating  oil,  fuel, 
hydraulic  fluid,  etc.  To  see  if  such  substances  could  cause  thermovideo- 
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graph  anomalies  that  might  ba  Interpreted  aa  flaws,  clean  lubricating  oil 
and  aollad  machine  grease  were  applied  to  the  abraded  Gr-D  sample  on  2.5 
cm  square  areas.  At  room  temperature  (Figure  9a),  a  clean  abrasion 
appears  about  0.2°C  cooler  than  unabradad  glossy  surface  as  usual,  while 
dean  lubricating  oil  does  not  change  surface  emlssivlty  characteristics 
from  thosa  of  the  glossy  paint.  In  fact,  the  oil  appears  to  "cover  up" 
the  slight  apparent  temperature  effect  of  the  abrasion.  Heating  for  one 
minute  with  a  radiant  heat  source  and  cooling  for  three  minutes  yields  a 
thermovideograph  where  no  apparent  difference  can  be  seen  between  glossy 
paint,  abrasion,  and  abrasion  coated  with  clean  lubricating  oil  (Figure  9b). 
Continued  cooling  for  a  total  of  6.5  minutes  gives  different  results 
(Figure  9c):  both  the  clean  abrasion  and  the  oil  coated  abrasion  appear 
0.1  to  0.2°C  cooler  than  surrounding  glossy  paint.  These  differences 
are  extremely  small,  however,  and  could  even  have  been  caused  by  nonuniform 
heating  or  cooling  of  the  sample. 

b.  A  photograph  of  sample  Gr-D  containing  a  clean  abrasion,  an 
abrasion  with  a  heavy  coat  of  soiled  machine  grease,  and  a  light  2.5  cm 
x  2.5  cm  coat  of  the  same  soiled  grease  on  virgin  painted  surface  is 
shown  in  Figure  10.  Thermovideographs  are  shown  in  Figure  11.  At  room 
temperature  (Figure  11a)  there  is  little  difference  from  the  surrounding 
area.  When  heated  one  minute  and  cooled  for  30  seconds, however,  the  heavy 
coat  of  soiled  grease  shows  a  minimum  of  2°C  apparent  temperature  difference, 
and  the  light  coat  of  soiled  grease  appears  nearly  1.4°C  higher  in  temperature 
than  the  clean,  painted  surface  (Figure  lib) . 

5.  DISCUSSION.  Surfaces  coated  with  glossy  Navy  aviation  paint  can 
appear  thermographlcally  different  than  virgin  aluminum  and  epoxy  matrix 
composites.  Paint  thickness  variations  are  clearly  delineated  on  aluminum 
surfaces,  but  are  not  thermographlcally  visible  on  epoxy  matrix  surfaces. 

The  high  gloss  of  the  Navy  aviation  paint  reflects  concentrated  heat 
sources  such  as  light  bulbs,  and  surfaces  to  be  inspected  must  be  shielded 
from  such  heat  sources.  Heavy  clean  abrasions  on  painted  surfaces  can 
appear  a  few  tenths  of  a  degree  centigrade  cooler  than  virgin  painted 
areas,  but  should  not  be  mistaken  for  flaws  because  the  difference  is  so 
slight.  Clean  lubricating  oil  does  not  appear  differently  from  clean. 
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painted  surfaces,  and  in  fact  seems  to  reduce  the  slight  effect  of  abrasions. 
Soiled  grease,  on  the  other  hand,  can  appear  several  degrees  centigrade 
higher  than  clean,  painted  surfaces. 

0.  THERMAL  FIELD  APPLICATION  METHODS 

1.  HEAT  GUN.  As  a  result  of  the  preliminary  studies  (reference  (a)), 
it  became  apparent  that  a  heat  gun  method  of  applying  energy  to  the  sample 
would  not  be  effective.  For  the  inspection  of  large  surface  areas,  the 
heat  gun  was  too  small  a  source.  In  addition,  this  method  provided 
limited  heat  input  into  the  sample.  The  following  other  methods  were 
examined: 

-  Heat  Lamp 

-  Local  application  of  strip  heaters 

-  Radiative  thermal  source 

2.  HEAT  LAMP.  Heat  lamps  have  the  potential  to  be  relatively  in¬ 
expensive  and  easy  to  maintain  in  the  naval  aviation  maintenance  environ¬ 
ment.  Ideally,  banks  of  heat  lamps  could  comprise  a  hand-held  heater; 
burn-outs  could  be  easily  replaced  from  a  stock  of  spares.  For  this 
study,  one  250-watt  heat  lamp  was  used  to  heat  selected  samples  of  graphite/ 
epoxy.  The  resulting  heat  pattern  was  circular  since  there  was  only  one 
bulb.  The  heat  lamp,  although  providing  adequate  heat  for  certain  NDE 
purposes,  was  judged  not  to  have  high  enough  heat  output  for  detection  of 
flaws  existing  deep  within  a  structure.  Perhaps  this  might  be  overcome  by 
using  higher  wattage  bulbs.  The  largest  drawback  to  the  use  of  the  heat 
lamp  is  the  tendency  of  examined  surfaces  to  reflect  the  heat  from  the 
bulb.  If  heat  lamps  are  to  be  used,  they  must  be  carefully  shielded  from 
surfaces  under  examination  to  avoid  unwanted  spurious  thermographic  out¬ 
put. 

3.  LOCAL  STRIP  HEATERS 

a.  The  use  of  local  surface  heaters  was  considered  because  of  the 
simplicity  and  economy  of  the  method.  A  technician  could  apply  the  heaters 
over  the  surface  to  be  inspected  end  then  examine  the  temperature  profiles 
with  the  thermovision  camera  as  the  energy  was  conducted  through  the  sample. 
The  results  could  then  be  compared  with  known  unflawed  standards.  The 
technique  would  be  based  upon  steady-state  conditions  and  require  only  one 
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technician.  (A  radiative  source  method  would  examine  transient  response 
and  therefore  would  require  two  technicians,  one  applying  the  thermal 
field  and  the  second  operating  the  camera.  The  advantages  of  this  method 
are  no  surface  contact  is  required,  and  the  thermal  field  application 
would  be  uniform.) 

b.  To  examine  the  effectiveness  of  the  local  field  application 
method,  a  15- cm  x  13- cm  x  . 24-cm  sample  of  graphite/epoxy  was  used. 

Figure  12  presents  a  sketch  of  the  sample  which  included  a  0.65-cm  diameter 
partial  through-hole.  The  sample  was  8-ply  unidirectional  laminate  and 
was  not,  therefore,  isotropic.  The  sample  was  mounted  in  a  manner  similar 
to  the  surface  finish  samples.  Heat  was  applied  by  a  heat  gun  in  one  case 
and  by  Minco  surface  heaters  mounted  directly  to  the  surface.  For  the 
test  using  the  heat  gun,  the  sample  was  heated  uniformly  and  then  allowed 
to  cool  by  free  convection.  In  the  case  of  the  surface  heater,  power  was 
applied  to  the  heater  to  establish  a  given  steady-state  condition.  The 
energy  was  conducted  along  the  sample  and  then  was  removed  by  free  convec¬ 
tion.  A  simplified  analytical  model,  which  considered  the  sample  to  be  a 
heated  fin,  was  also  examined  so  that  comparisons  with  the  experiments 
could  be  made. 

c.  Analytical  prediction  of  the  thermal  response  of  the  graphite/ 
epoxy  sample  with  the  heater  applied  to  the  surface  was  obtained  using  a 
quasi  one-dimensional  fin  approximation.  In  this  approach,  the  sample  is 
assumed  to  have  temperature  gradients  in  the  vertical  direction  only.  The 
sample  was  then  divided  into  two  regions;  a  heated  region  and  an  unheated 
region.  The  governing  equations  for  each  region  are  given  below  (reference 
<b»: 

Region  A  (Heated) 

dx 

where: 

Ta  “  Temperature  in  Region  A 

Q  -  Heat  dissipated  by  heater  per  unit  cross  sectional  area 
h  "  Convective  heat  transfer  coefficient 
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P  ■  Perimeter  of  fin 
A  ■  Cross  sectional  area  of  fin 
Tw  ■  Ambient  temperature 
k  "  Thermal  conductivity 
x  "  Position  along  sample 


Region  B  (Unheated) 


B  hP 


<tb-  T«) 


(I 1-2) 


where: 


*  Temperature  in  Region  B 


The  boundary  conditions  are: 
x  -  0:  dTA/dx  -  0 
x  -  a:  T.  -  T-,  dT. 


b:  dT_ 


Solving  these  equations  simultaneously  and  applying  the  boundary  conditions, 
equations  (11-3)  and  (II-4)  are  obtained. 


^B  T<®  +  [a inh  ma  cosh  mx  -  sinh  mx] 

tanh  mb 

T.  ■  T  +  <JA  r,  .  sinh  m(a-b)  cosh  mxq 
hP  1  sinh  nb  3 


(1 1—3) 


(II-4) 


0  <  x  <  a 

where:  m  »  (hP/kA)1* 

Using  equations  (II -3)  and  (II-4),  the  temperature  distribution  in  the  fin 
(sample)  as  a  function  of  position  for  different  values  of  h,  k,  and  Q  can 
be  obtained. 

d.  Figure  13  presents  the  results  for  the  case  of  the  graphite/ 
epoxy  sample  exposed  to  an  external  thermal  field  applied  by  a  heat  gun.  In 
Figure  13a,  the  heat  input  was  applied  vertically,  and  in  Figure  13b,  the 
heat  input  was  horizontal.  It  should  be  noted  that  the  thermal  conductivity 
was  much  greater  in  the  horizontal  direction  than  in  the  vertical  direction. 
The  white  lines  in  the  figures  are  particular  Isotherms  which  the  camera 
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identifies.  In  both  Figure*  13*  and  13b,  the  isotherm  is  distorted  in  the 
middle  of  the  picture  where  the  non- through  hole  defect  is  located.  In 
addition,  the  thermal  pattern  is  distorted  at  the  top  of  the  picture  where 
an  edge  delamination  has  been  introduced  into  the  Bample.  Figure  13c  pre¬ 
sents  resultB  for  the  same  case  as  13a  and  13b  but  with  sensitivity  set  at 
20°C  from  black  to  white  rather  than  50°C.  The  center  defect  is  more 
pronounced  because  of  the  small  range  between  shades. 

e.  The  results  for  the  case  in  which  the  heat  was  applied  by  means 
of  a  surface  heater  are  entirely  different.  Figures  14a,  14b,  and  14c  pre¬ 
sent  results  for  the  cases  of  heat  inputs  of  1,  4,  and  10  watts  respectively, 
Heat  flows  from  beneath  the  heater  by  conduction,  and  the  flow  pattern  is 
altered  when  a  defect  is  encountered.  For  each  case  in  Figure  14,  the 
temperature  gradient  is  so  steep  that  no  discernible  temperature  effects 
reach  the  defect  in  the  sample.  In  Figure  14a,  the  nonuniformities  near 

the  top  of  the  picture  are  due  to  heater  nonuniformities.  The  steep 
gradients  are  due  to  the  low  thermal  conductivity  in  the  vertical  direction. 
Use  of  this  approach  for  anisotropic  laminates  would  require  heater  applica¬ 
tion  near  the  defect.  Since  in  the  HOT  situations  the  flaw  location  is 
not  known,  this  approach  would  not  be  acceptable. 

f.  Theoretical  results  for  the  anisotropic  case  are  shown  in  Figure 
IS.  They  agree  with  the  experimental  evidence  that  heat  flow  in  the 
vertical  direction  is  negligible.  In  Figure  14,  the  sensitivity  is  5°C, 
meaning  that  the  temperature  range  from  black  to  white  is  5°C.  During  the 
test,  the  black  level  was  maintained  at  room  temperature  so  that  the  middle 
of  the  picture  in  Figure  14a  is  about  1°C  above  room  temperature.  This  is 
the  approximate  value  at  X/L  -  0.5  and  a  heat  input  of  1  watt  in  Figure  15. 
Similar  agreement  exists  at  10  watts.  Unfortunately,  a  quasi-isotropic 
laminate  of  graphite/ epoxy  was  not  available  at  the  time  of  these  tests. 
However,  a  computer  run  for  an  isotropic  sample  indicated  that  the  gradient 
is  reduced  considerably. 

4.  RADIANT  HEAT  SOURCE 

a.  Quasi-isotropic  [0/±45/90]g  graphite/epoxy  samples  were  obtained 
for  the  radiative  field  application  evaluation.  Initially,  the  1320-watt 
space  heater  shown  in  Figure  16  was  used.  The  samples  used  were  8-ply  and 
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32-ply  quasi-isotropic  graphite/ epoxy.  Each  sample  was  30.8  cm  x  46  cm 
and  contained  six  manufactured  delaminations.  The  delaminations  were  6.4 
mm  x  6.4  mm,  12.7  nm  x  12.7  mm,  and  25.4  mm  x  25.4  mm  and  were  located  at 
two  different  depths  through  the  sample.  This  arrangement  allowed  for  both 
size  and  depth  detection  studies.  The  samples  were  mounted  vertically  as 
shown  in  Figure  17. 

b.  Initial  studies  with  the  space  heater  and  8-ply  sample  were 
very  encouraging.  This  method  was  preferable  to  the  heat  gun  method  which 
has  some  bias  because  defect  location  was  known.  However,  with  the  thicker 
32-ply  sample,  the  commercial  space  heater  proved  to  be  an  inadequate  source. 
For  the  8-ply  sample,  the  defect  was  clearly  visible  for  2  minutes.  At 

32  plies,  there  was  disagreement  between  experimenters  if  the  defect  was 
visible  at  all.  Those  who  contended  that  the  defect  was  visible,  however , 
agreed  that  it  was  only  for  1  or  2  seconds.  For  field  use,  this  was 
minimally  satisfactory.  These  results  indicated  that  the  heat  source  was 
not  large  enough  to  cause  significant  thermal  effect  at  the  defect.  Quasi¬ 
isotropic  graphite/ epoxy  has  a  large  lateral  conductivity  (40  W/m°C)  com¬ 
pared  to  the  thermal  conductivity  normal  to  the  sample  (0.5  W/m°C) .  To 
cause  significant  effects  through  the  sample  thickness,  a  large  thermal  flux 
normal  to  the  surface  must  be  provided. 

c.  To  determine  the  required  heat  input  to  the  sample,  a  numerical 
model  of  the  sample  was  created  using  50  nodes,  and  was  analyzed  by  the 
heat  transfer  code  described  earlier.  The  physical  size  of  the  model  was 
12.7  cm  x  12.7  cm.  There  were  25  surface  nodes  exposed  to  the  flux  and 

25  rear  nodes.  The  conductive  couplings  between  nodes  accounted  for  the 

difference  in  normal  and  lateral  thermal  conductivity  of  graphite/ epoxy . 

A  delamination  was  simulated  by  varying  the  coupling  between  the  center 

surface  node  and  the  rear  node.  Figure  18  presents  results  for  the  32- 

ply  graphite/epoxy  laminate  with  25.4-mm  and  12.7-mm  delaminations  at 

the  center  node.  The  results  show  temperature  of  the  sample  3  minutes  after 

2 

exposure  to  a  0.363  W/cm  flux.  Depending  upon  the  size  of  the  delamination, 
a  temperature  distortion  of  0.25  to  1.0°C  will  exist  for  this  flux.  This 
magnitude  is  considered  within  the  capability  of  the  thermovision  camera. 

A  flux  of  0.363  W/cm  is  very  large  and  is  equivalent  to  approximately  3 
earth  solar  constants.  Using  this  data,  a  new  radiative  heat  source  was 
designed. 
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d.  Initially,  the  design  considered  a  line  source  located  a  fixed 
distance  from  the  surface.  Figure  19  presents  the  geometry  examined  in  the 
design.  The  parameters  in  the  design  were  the  diameter  of  the  source,  dis¬ 
tance  from  the  surface,  and  temperature  of  the  source.  The  goal  of  the 
design  was  to  maintain  the  required  flux  at  a  distance  of  15  cm  of  the 
source  centerline.  Figure  20  presents  a  comparison  of  source  temperature 

versus  source  diameter  for  different  distances  from  the  surface  at  a  flux 
2 

of  0.363  W/cm  .  From  this  figure,  it  is  apparent  that  a  single  source 

requires  excessive  temperatures  for  a  reasonable  size.  An  alternate 

approach  is  to  use  multiple  sources  each  providing  a  lower  flux.  For 

2 

Instance,  two  sources  providing  0.187  W/cm  each  would  operate  at  a 
significantly  lower  temperature  and  provide  the  same  input  to  the  surface. 
The  final  design  was  based  upon  a  two-source  configuration.  Figure  21 
presents  the  geometry  of  the  final  radiative  heat  source  design.  Figure 
22  is  a  photograph  of  the  actual  heater.  The  heating  elements  consist 
of  1.25-cm  glass  tubes  around  which  nichrome  wire  is  wound.  The  wire 
spacing  is  adjusted  to  maintain  a  uniform  temperature  field  on  the  rod. 

The  sides  of  the  heater  box  are  lined  with  polished  aluminum  to  reflect 
more  radiant  energy  onto  the  sample  surface.  The  use  of  reflected  energy 
reduces  the  required  power  of  the  heater  and  wire  temperatures. 

e.  The  high  powered  radiative  source  was  used  on  the  32-ply  samples 
with  excellent  results.  These  results  are  presented  in  the  next  section. 

E.  EVALUATION  OF  FLAW  DETECTION  CAPABILITIES.  It  has  already  been 
established  that  thermography  can  detect  flaws  in  composite  materials 
(reference  (a)).  The  purpose  of  this  section  is  to  describe  the  size, 
range,  and  types  of  flaws  which  can  be  located  with  the  AGA  thermovision 
system. 

1.  GRAPHITE/EPOXY  LAMINATES.  Graphite/epoxy  was  used  for  the  first 
portion  of  the  flaw  detection  studies.  Two  thicknesses  were  considered: 

8  ply,  1.19  mm.;  and  32  ply,  5.26  mm.  The  laminates  tested  were  a  quasi¬ 
isotropic  [0/±45/90]g  composite  with  approximately  an  80:1  ratio  between 
the  in-plane  and  normal- to- the-plane  thermal  conductivities;  and  a 
[02/±45]b  anisotropic  layup  with  a  4:1  ratio  of  major  to  minor  in-plane 
conductivities;  and  approximately  160:1  ratio  between  major  in-plane  and 
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normal- to- the- plane  conductivities.  Two  types  of  flews  were  examined, 

|  namely  simulated  delaminations  and  non- through  holes. 

a.  8-Ply  Graphite/Bpoxy  Laminates 

(1)  The  simulated  delaminations  consisted  of  rectangular  bags 

*  of  hollow  glass  beads  introduced  into  the  sample  during  manufacturing. 
These  square  delaminations  varied  In  size  from  6.35  mm  to  25  mm  square. 

i  The  flaws  were  located  1/8  of  the  thickness  into  the  sample  and  1/2  of 

the  thickness.  Figure  17  shows  the  8-ply  sample  mounted  in  the  test 
fixture. 

(2)  The  non-through  hole  defects  which  simulate  subsurface 
shattering  due  to  impact  damage  were  introduced  into  undamaged  samples 

by  drilling  various  size  holes.  Both  the  diameter  and  depth  of  the  holes 
were  varied.  Hole  diameters  ranged  from  12.7  mm  to  6.4  am.  These  tests 
were  conducted  only  in  the  32-ply  samples  and  hole  depth  did  not  exceed 
one  half  of  the  sample  thickness. 

(3)  The  samples  were  mounted  vertically  as  shown  in  Figure 
23.  Figure  24  shows  the  arrangement  of  the  camera  and  the  labelling  of 
the  test  panels.  Because  of  optical  limitations,  the  camera's  field  of 
view  could  not  cover  the  entire  sample.  The  tests  were,  therefore, 
arranged  to  examine  two  flaws  at  a  time. 

(4)  Figure  25  presents  results  for  the  8-ply  sample  with  the 
external  field  applied  by  the  space  heater.  The  delaminations  are  clearly 
visible  in  the  photograph.  The  ten-color  scale  along  the  bottom  of  the 
picture  shows  the  temperature  range  of  the  sample.  This  scale  can  be 
adjusted  from  1°C  black  to  white  to  1000°C  black  to  white.  The  most 
effective  range  used  in  the  present  study  was  between  2  and  20°C. 

b.  32-Plv  Graphite/Epoxy  Laminates 

(1)  As  discussed  earlier,  the  space  heater  was  ineffective  on 
the  32-ply  sample.  The  heater  shown  in  Figure  22  was  used  on  the  32-ply 
samples.  For  all  the  tests,  the  heater  was  operated  at  1.0  KW  and  exposed 
to  the  sample  for  60  seconds.  The  Infrared  camera  lens  was  positioned 

|  approximately  2  meters  from  the  sample. 

i 

(2)  Figure  26  presents  results  showing  flaws  B2  and  C2  approx- 

i  imately  60  seconds  after  the  heat  source  was  removed.  The  camera  sensi- 

•  tivlty  level  is  10°C  black  to  white.  It  should  be  noted  that  for  this 
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32-ply  sample  che  defects  are  not  visible  to  the  naked  eye.  At  the  moment 
of  removal  of  the  heater,  the  camera  indicated  that  the  sample  had  been 
uniformly  heated.  The  sample  then  cools  by  free  convection  uniformly 
except  in  the  vicinity  of  the  defects.  Figure  27  presents  the  same  results 
for  the  case  in  which  flaws  A2  and  B2  are  Illuminated.  Flaw  A2  is  the  6.4 
mm  square  located  half  way  through  the  sample.  It  is  located  at  the  large 
white  spot  in  Figure  27a.  Flaw  B2  is  the  smaller  white  spot  to  the  right 
of  the  picture.  Figure  27b  was  taken  sometime  later  in  the  cooling  cycle. 
The  white  line  outlining  the  colored  region  is  an  isotherm.  Flaw  A2  is 
located  at  the  light  orange  spot  in  the  center  of  Figure  27b. 

(3)  Figure  28  presents  a  sequence  of  pictures  with  the  sample 
inverted  so  that  the  flaws,  located  1/8  of  the  way  through  the  sample,  are 
in  the  field  of  view  of  the  camera.  In  Figure  28a,  flaw  B1  is  clearly 
visible  at  the  right  of  the  picture.  Flaw  Cl  is  masked  in  the  hot  area 

at  the  center  of  the  picture.  It  is  difficult  to  get  both  flaws  (Cl  and  Bl) 
at  the  same  temperature  so  that  they  might  be  visible  early  in  the  cooling 
process.  Views  b,  c,  and  d  of  Figure  28  present  other  attempts  to  capture 
both  flaws  by  waiting  longer  in  the  cooling  period.  Figure  28c  represents 
the  best  results  clearly  showing  both  flaws  although  Cl  is  still  hotter. 

(4)  Views  a,  b,  and  c  of  Figure  29  present  a  sequence  of  pic¬ 
tures  of  flaws  A2,  B2,  and  C2  which,  because  they  are  at  the  midplane  of 
the  sample,  have  the  same  relative  position  to  the  surface  as  Figures  26 
and  27.  Figure  29a  shows  defect  C2  at  the  center  of  the  picture.  Figure 
29b  shows  C2  as  a  bright  spot  at  the  center  and  B2  at  the  right.  Figure 
29c  presents  the  smallest  defect  A2  at  the  right  side  of  the  picture. 

(5)  Following  these  tests,  the  sample  was  reversed  so  that 
flaws  Al,  Bl,  and  Cl  were  7/6  of  the  sample  thickness  from  the  heated 
surface.  This  represents  the  most  severe  test  of  the  flaw  detection 
capable  of  the  sample.  Figure  30a  and  b  shows  flaws  Al  and  Bl  clearly. 

The  camera  sensitivity  was  set  at  10°C  for  these  pictures.  Figure  30b 
shows  the  smallest  defect  Al  to  the  left  outlined  by  an  isotherm.  Figure 
31  presents,  for  the  same  sample  orientation,  defects  Bl  and  Cl.  Flaw  Cl 
appears  to  the  right  side.  Figure  31a  is  for  a  sensitivity  of  10°C  while 
for  Figure  31b,  the  setting  was  5°C.  These  results  show  that  if  a  higher 
power  heat  source  is  used,  flaws  in  thicker  samples  can  be  located.  Un¬ 
fortunately,  the  32  ply  was  the  thickest  specimen  available. 
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(6)  All  of  the  preceding  results  were  obtained  on  the  quasi- 
isotropic  [Q/±45/90jg  graphite/epoxy  samples.  In  order  to  test  the  effects 
of  in-plane  thermal  anisotropy,  partial  through-holes  were  drilled  in  an 
anisotropic  graphite/epoxy  laminate  composed  of  32  plies,  and 

having  in-plane  principal  thermal  conductivities  of  80W/m  deg  C  and  20  W/m 
deg  C.  Hole  diameters  were  small  (6  to  11  mm)  and  ranged  from  25Z  to  87Z 
through  the  thickness.  The  anisotropic  specimen  was  heated  and  observed 
from  the  intact  surface;  no  evidence  of  the  partial  through-holes  could 
be  seen  with  the  unaided  eye.  Figure  32  presents  the  defect  map  for  the 
specimen  as  seen  from  the  intact  side.  Mote  the  existence  of  an  un¬ 
intentional  delamination  in  the  sample,  at  an  unknown  depth,  which  was 
identified  by  the  manufacturer's  C-scan. 

C7)  Defects  were  examined  two  or  three  at  a  time  since  the 
Thermovision  could  not  accommodate  the  entire  sample  in  its  field  of  view. 
Figure  33  presents  typical  thermovideograph  results  after  heating  with  a 
radiant  heat  source  for  10  to  30  seconds,  and  cooling  for  5  to  10  seconds. 

In  Figure  33a,  a  6-mm-dla,  87Z  through-hole  shows  as  a  large  bright  spot 
to  the  left,  and  an  11-mm-dia,  65Z  through-hole  can  be  seen  to  the  right. 

An  11-mm-dia,  25Z  through-hole  centrally  located  above  the  two  bright 
spots,  at  the  apex  of  an  equilateral  triangle,  cannot  be  observed.  In 
fact,  neither  of  the  25Z  through-holes  could  be  identified  in  this  sample. 
Figure  33b  shows  an  11-mm-dia,  65Z  through-hole  to  the  left,  and  an  11-mm- 
dia,  50Z  through-hole  to  the  right,  at  an  infrared  camera  sensitivity  of 
2°C.  The  defect  nearer  the  observed  surface  appears  more  clearly,  as 
expected.  Mote  that  the  infrared  camera  also  shows  a  hot  spot  between  the 
two  man-made  flaws,  corresponding  to  the  size  and  location  of  the  afore¬ 
mentioned  unintentional  delamination.  Figure  33c  shows  the  same  three  flaws 
at  a  5°C  camera  setting.  Finally,  it  is  noted  that  the  laminate's  thermal 
anisotropy  did  not  appear  to  affect  the  detectability  of  the  defects;  the 
only  difference  between  the  an<*  [0/±45/90]a  results  is  a  slight 

elongation  of  the  "hot  spots"  in  the  direction  of  greatest  conductivity 
for  the  anisotropic  sample. 


2.  S2-GLASS/EP0XY  LAMINATES 

a.  Glass/epoxy  composites  are  common  in  several  different  types  of 
aircraft.  Since  the  in-plane  thermal  conductivities  of  glass/epoxies  are 
usually  an  order  of  magnitude  or  two  less  than  graphite/epoxies  (about 

0.5  W/m°C  versus  40  W/m°C  for  quasi-isotropic  laminates),  it  was  judged 
necessary  to  test  a  typical  glass/epoxy  material  for  flaw  detection  ability. 
An  S2-glas8/epoxy,  3M  SP-250-S2,  was  chosen.  Panels  up  to  16  ply  having  a 
thickness  of  about  4  mm  were  fabricated  with  plastic  bags  (empty)  to  simu¬ 
late  flaws.  Unfortunately,  the  bags  bonded  during  the  curing  process  and 
did  not  provide  a  thermal  image  when  heated.  Therefore,  the  laminates 
were  Impacted  by  a  blunt  object  several  times  at  a  velocity  of  about  3  m/sec 
to  open  the  delaminations.  Only  the  2.5  cm  square  delaminations  were  able 
to  be  opened  in  the  32-ply  sample.  In  addition  to  the  simulated  delamina¬ 
tions,  partial- through-holes  from  6  mm  to  13  mm  in  diameter  and  up  to  802 
through  the  thickness  were  drilled  in  the  samples.  Figure  34  is  a  defect 
map  for  the  16-ply  [0/±45/90]  panel. 

b.  Results  of  testing  the  16-ply  panel  are  shown  in  Figure  35. 

Two  partial-through-holes,  one  20  mm  dia,  152  through  (left)  and  one  11  mm 
dla,  602  through  (right)  are  seen  above  a  25-nm-square  delamination  one 
ply  from  the  top,  10  seconds  after  heating,  in  Figure  35a.  Figure  35b 
shows  the  same  section  of  panel  after  cooling  for  six  minutes.  Nbte  that 
the  free  convection  cooling  has  reversed  the  temperature  of  the  surface 
over  the  holes,  and  that  the  delamination  is  no  longer  visible.  A  25  mm 
square  delamination  located  halfway  through  the  thickness  was  heated  for 
10  seconds  and  cooled  for  10  seconds,  and  resulted  in  the  thermovideograph 
of  Figure  35c.  The  delamination  is  visible  by  its  0.3°C  temperature  rise 
above  surrounding  area.  In  Figure  35d,  three  partial-through-holes  are 
visible  as  hot  spots  (lef t-to-rlght:  20  mm  dia,  152  through;  11  mm  dia, 

602  through;  and  11  mm  dia,  302  through). 

c.  Four  other  delamination-type  defects  could  not  be  detected; 
these  were  two  (0.5  cm  square)  implanted  one  ply  and  8  plies  down  from 
the  surface,  respectively,  and  two  others  induced  by  repeated  impact  with 

a  blunt  object.  The  two  implanted  delaminations  were  not  true  separations, 
as  previously  indicated.  The  type  and  extent  of  the  blunt  object-induced 
damage  sites  is  not  known.  Therefore,  it  is  impossible  at  this  time  to 
draw  meaningful  conclusions  from  these  four  defects. 
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3.  DISCUSSION .  The  flaw  detection  capability  studies  have  shorn 
that  delaminations  and  partial-through  holes  can  be  detected  in  quasl- 
isotropic  and  anisotropic  laminates  of  graphite/epoxy  and  S2-glass/ epoxy. 
Sizes  ranged  from  0.6  cm  to  2.5  cm  across,  located  as  far  as  7/8  through 
the  thickness  of  up  to  32  plies  thick  graphite/epoxy  and  16  plies  thick 
S2-glass/ epoxy.  The  infrared  camera  was  even  able  to  detect  an  uninten¬ 
tional  flaw  induced  during  manufacturing.  Anisotropy  does  not  seem  to  have 
any  significant  effect  on  the  detection  of  flaws.  The  thickest  laminates 
available  were  5-mm  32-ply  graphite/epoxy  and  4-zzm  16-ply  S2-glass/ epoxy 
samples;  it  does  not  appear  that  the  limits  of  size  and  location  on 
delamination  detection  capability  have  yet  been  found. 
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III.  STRESS-GENERATED  THERMAL  FIELD  METHOD 

A.  OBJECTIVES  AND  OUTLINE  OF  THE  SGTF  STUDY 

1.  The  SGTF  study  had  as  Its  primary  objective  the  determination  of 
the  feasibility  of  detecting  flaws  in  composite  and  metal  structures  by 
generating  heat  near  the  flaws  with  high-frequency,  low-stress  cyclic 
loading.  It  had  been  known  that  partial-through  holes,  notches,  and 
delaminations  can  be  detected  by  SGTF  methods  in  glass/epoxy,  boron/ 
epoxy,  and  graphite /epoxy  structures.  It  is  not  known,  however,  whether 
composite  structures  can  be  damaged  by  the  levels  of  cyclic  stress  applica¬ 
tion  necessary  to  generate  heat  for  thermal  flaw  detection.  Also,  SGTF 
techniques  have  not  been  successful  for  metal  structures.  The  planned 
SGTF  study  was  to  first  determine  the  levels  of  stress  and  frequency  which 
can  be  applied  to  structural  composites  and  metals  without  causing  struc¬ 
tural  damage  and  still  obtain  a  thermal  image  allowing  flaw  detection. 

2.  The  study  begins  by  first  performing  a  semi -empirical  analysis  of 
the  heat  generated  in  cyclically  loaded  composites,  and  experimentally 
calibrating  two  materials — a  graphite /epoxy  and  an  S2-glasa/ epoxy  com¬ 
posite.  The  results  of  the  semi -empirical  heat  generation  analysis  are 
then  used  in  conjunction  with  a  finite  element  stress  analysis  to  predict 
heat  generated  near  a  6-nmr  dia  hole  in  graphite /epoxy  and  glass/epoxy 
tensile  specimens.  The  analytical  results  are  compared  to  the  results 

of  tests  run  on  these  specimens.  Aluminum  specimens  were  also  tested  and 
their  results  are  also  described.  Finally,  a  feasibility  assessment  of 
SGTF  techniques  in  composites  and  metals  is  made,  and  the  results  of  a 
preliminary  investigation  into  ultrasonic  SGTF  methods  is  presented. 

B.  ANALYSIS  OF  HEAT  GENERATION  IN  CYCLICALLY  LOADED  COMPOSITES 

1.  APPROACH  TO  THE  PROBLEM 

a.  Determination  of  heat  generation  in  composite  materials  under 
cyclic  loading  can  be  approached  in  several  ways.  At  the  one  extreme,  a 
detailed  analysis  of  viscoelastic  stresses  in  a  unidirectional  composite 
could  be  made,  and  the  heat  energy  per  unit  time  dissipated  due  to  visco¬ 
elastic  behavior  could  be  calculated.  This  would  necessitate  knowing  the 
in  situ  viscoelastic  properties  of  the  epoxy  matrix  material.  Even  if 
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an  exCensive  testing  program  could  be  undertaken  to  determine  these  prop- 
artles,  it  is  likely  that  the  average  in  situ  matrix  properties  would  not 
be  the  same  due  to  curing  differences  and  the  existence  of  voids,  imper¬ 
fections,  disbonds  between  fiber  and  matrix,  etc.  At  the  other  extreme 
lies  a  purely  empirical  approach  which,  through  cyclic  thermomechanical 
testing,  could  experimentally  determine  the  viscoelastic  heat  generated 
for  a  wide  range  of  composite  laminates,  combined  stress  states,  and 
cyclic  frequencies.  This  totally  empirical  approach  is  unrealistic  be¬ 
cause  of  the  vast  number  of  specimens  and  tests  which  would  be  required. 

b.  The  approach  adopted  for  this  study  chose  appropriate  elements 
from  both  extremes  combined  into  a  semi-empirical  analysis .  First,  elastic 
stresses  were  determined  in  a  representative  volume  element  of  unidirectional 
composite  under  a  general  state  of  in-plane  combined  stress.  Next,  a 
suitable  viscoelastic  material  model  was  chosen,  and  the  form  of  an 
equation  for  viscoelastic  point  heat  energy  generated  was  developed  from 
this  model.  It  was  assumed  that  the  elastic  and  viscoelastic  stress  states 
were  not  significantly  different,  and  the  elastic  stress  distribution  was 
combined  with  the  viscoelastic  point  heat  energy  model  to  give  an  equation 
for  average  heat  energy  dissipated  in  a  unidirectional  layer  of  composite. 

The  equation  contains  two  constants  which  can  be  determined  for  a  given 
composite  material  system  from  a  limited  number  of  tests  on  simple  laminates. 

c.  With  this  semi-empirical  equation  for  heat  generated  in  a  uni¬ 
directional  composite,  the  procedure  for  determining  heat  generated  and 
temperatures  in  a  laminated  structural  composite  is  as  follows:  First, 
a  global  structural  elastic  stress  analysis  is  performed  under  the  cyclic 
loads  using  any  conventional  method  (such  as  finite  element  analysis). 

Next,  the  structure  is  divided  into  convenient  subregions,  and  stresses  in 
each  layer  are  calculated  using  laminate  analysis  techniques.  The  semi- 
emplrlcal  heat  generation  equation  is  then  used  for  each  layer  of  each 
subregion  to  determine  rate  of  heat  generation.  Finally,  a  heat  transfer 
analysis  is  performed  on  the  discretized  structure,  and  temperature  distri¬ 
butions  are  calculated.  Figure  36  illustrates  this  procedure. 
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2.  STBMSgg  IN  UNIDIRECTIONAL  COMPOSITES 

e.  In  order  to  develop  e  semi-empirical  equation  for  beet  generated 
under  cyclic  loading  in  composites,  it  is  desirable  to  know  the  detailed 
state  of  stress  and  strain  in  the  fiber  and  matrix  phases  due  to  cyclically 
applied  unidirectional  layer  average  stresses.  The  heat  energy  generated 
due  to  viscoelastic  behavior  can  then  be  calculated  by  Integrating  stress 
times  differential  strain  over  one  cycle  of  applied  layer  average  stresses 
and  dividing  by  the  time  for  one  load  cycle.  The  exact  determination  of 
viscoelastic  stresses  and  strains  in  the  fiber  and  matrix  phases  of  a 
composite  layer  is  a  very  difficult  task,  and  requires  knowledge  of  fiber 
and  matrix  viscoelastic  properties.  Nearly  all  high  strength  fibers  can 
be  considered  elastic  compared  to  epoxy  matrix  materials.  Hence,  nearly 
all  heat  generated  in  epoxy  matrix  composites  will  be  due  to  the  visco¬ 
elasticity  of  the  matrix  material.  Therefore,  it  is  assusied  herein  that 
fibers  are  purely  elastic,  and  that  all  heat  generated  will  be  in  the 
matrix.  Even  with  this  assumption,  however,  exact  solution  to  the  visco¬ 
elastic  stresses  and  strains  is  very  difficult.  In  composites.  In  the  part 
of  the  stress-strain  range  which  is  below  the  point  where  damage  to  either 
fiber  or  matrix  can  occur,  the  matrix  is  elastic  enough  to  calculate  stresses 
for  engineering  purposes  using  a  linear  elastic  analysis.  It  will  there¬ 
fore  be  assumed  in  this  study  that  stresses  (not  strains)  in  the  epoxy 
matrix  of  a  unidirectional  layer  can  be  determined  in  terms  of  the  applied 
average  layer  stresses  by  an  elastic  analysis.  Resulting  strains  can  then 
be  determined  from  the  viscoelastic  constitutive  relations.  It  is  noted 
that  stress  equilibria  will  be  satisfied  by  this  procedure,  but  not 
necessarily  compatibility. 

b.  Determination  of  matrix  stresses  in  a  unidirectional  composite 
layer  requires  choice  of  fiber  and  matrix  phase  geometry.  Actual  distribu¬ 
tion  of  fibers  in  the  matrix  is  statistical  in  nature,  with  some  regions 
approximated  by  a  rectangular  array  and  others  approximated  by  an  (irregular) 
hexagonal  array  (Figure  37).  These  deterministic,  regular  fiber/matrix 
geometries,  however,  make  it  difficult  to  obtain  analytical  stress  distri¬ 
butions  by  other  than  numerical  analyses  (references  (c)  and  (d)).  The 
randomness  of  fiber  distribution  which  exists  in  a  composite  layer  can  be 


30 


NAEC-92-157 


approximated  by  the  Hashin-Rosen  composite  cylinder  assemblage  (CCA) 
(reference  (e))  which  assumes  that  the  composite  consists  of  composite 
cylinders;  each  composite  cylinder  consisting  of  an  inner  fiber  cylindri¬ 
cal  core  surrounded  by  a  matrix  annulus  (Figure  37c).  Each  composite 
cylinder  may  have  a  different  outer  diameter,  but  all  have  the  same  ratio 
of  fiber  radius  to  composite  cylinder  outer  radius,  a/b.  The  volume 
fraction  of  fibers  in  the  composite  is  therefore  always  spacewise  constant 
on  an  average  basis.  Use  of  the  CCA  geometry  approximation  has  several 
advantages:  It  has  been  demonstrated  that  average  elastic  composite  prop¬ 
erties  determined  using  CCA  geometry  are  very  close  to  experimental  results 
and  numerical  results  using  rectangular  or  hexagonal  arrays.  In  addition, 
closed  form  equations  for  displacements  have  been  determined  by  Hashln  and 
Rosen  (reference  (e)  and  (f))  from  which  matrix  stresses  can  readily  be 
determined  for  all  but  one  case  of  average  layer  stresses.  Therefore,  the 
CCA  phase  geometry  approximation  is  adopted  here. 

c.  Consider  a  unidirectional  layer  of  composite  material  under 
in-plane  loading  (Figure  38).  The  axial  (1)  direction  is  parallel  to  the 
fibers,  the  transverse  (2)  direction  is  perpendicular  to  fibers  in  the 
plane  of  the  layer,  and  the  normal  (3)  direction  is  perpendicular  to  the 
plane  of  the  layer.  Average  in-plane  composite  layer  stresses  are  axial 
normal  stress,  transverse  normal  stress,  a 22 5  and  axial  shear  stress, 

0^2"  The  matrix  material  is  assumed  to  be  isotropic,  and  the  fibers  are 
assumed  to  be  transversely  isotropic  with  principal  elastic  axis  along 
fiber  axis.  Hashln  (reference  (f))  gives  displacements  for  the  cases  of 
axial  normal  straining,  axial  shearing,  and  plane  strain  isotropic 
straining  in  the  2-3  plane.  He  also  gives  displacement  and  stress  fields 
for  calculating  upper  and  lower  bounds  to  composite  shearing  modulus  in 
the  2-3  plane,  but  these  latter  displacement  fields  are  not  exact  solutions. 
There  are  no  closed  form  exact  solutions  to  transverse  normal  straining  of 
a  composite. 

(1)  Axial  Normal  Stress.  Stresses  in  matrix  material  due  to 
axial  straining  of  a  unidirectional  composite  are  given  in  cylindrical 
polar  coordinates  (Figure  37d)  as: 
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Where  a  denotes  stress,  r  Is  radial  coordinate,  0  is  hoop  coordinate,  z  Is 
axial  coordinate,  superscript  and  subscript  m  denotes  matrix,  K  and  G  are 
bulk  and  shearing  modulus,  respectively,  e°  Is  applied  layer  axial  strain¬ 
ing,  L  and  N  are  constants  which  are  functions  of  matrix  elastic  con- 

ID  ID 

stants,  and  and  are  constants  of  Integration  which  are  found  from 
interface  and  boundary  conditions  (see  reference  (f)  for  more  detail).  Since 
the  composite  is  in  uniaxial  straining,  the  average  axial  normal  composite 
stress  is  related  to  the  applied  axial  strain  e°  by 
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where  is  the  effective  axial  modulus  of  the  unidirectional  layer. 
Using  (III-2)  and  the  boundary/continuity  conditions,  the  following 
equations  are  obtained  for  matrix  stresses  due  to  average  axial  normal 


composite  layer  stress  (reference  (g)): 
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with  remaining  stresses  equal  to  zero  (see  eq . (III-ld,e,f) ,  and 
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Here,  E  is  Young's  modulus,  v  is  Poisson's  ratio,  V  is  volume  fraction, 
subscript  £  refers  to  fiber  quantities,  and  remaining  symbols  are  as 
before. 

(2)  Axial  Shear.'  If  u^,  u^,  u^  represent  displacements  In  the 
layer  coordinate  system  (Figure  38) ,  then  applying  the  following  displace¬ 
ments  to  the  lateral  composite  cylinder  surface  is  equivalent  to  macro¬ 
scopic  shear  straining: 
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with  being  the  macroscopic  tensorial  shear  strain  (a  constant).  These 
displacements  give  rise  to  matrix  stresses  of  (see  reference  (g)): 
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where  the  composite  layer  average  shear  stress  given  by 
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is  the  effective  shearing  modulus  of  the  unidirectional  layer,  and 
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Y-Gf/GB.  (III-5j) 

(3)  Transverse  Normal  Stress.  As  previously  stated,  no  exact 
solution  exists  for  matrix  stresses  In  a  unidirectional  composite  layer 
under  transverse  normal  straining  or  stressing.  Hashin  (reference  (f)), 
however,  has  presented  the  displacement  field  solution  for  plane  isotropic 
straining.  He  also  presents  a  displacement  field  and  a  stress  field  used 
to  obtain  upper  and  lower  bounds,  respectively,  to  the  effective  composite 
transverse  shear  modulus.  While  neither  of  these  “upper  and  lower  bound" 
solutions  has  been  shown  to  be  an  exact  solution  to  the  problem  of  applied 
composite  transverse  shear  stress,  it  is  believed  that  the  displacement 
field  used  for  the  upper  bound  on  effective  transverse  shear  modulus  is 
closer  to  the  exact  solution.  Therefore,  the  “upper  bound"  displacement 
field  is  used  here  to  obtain  matrix  stresses  for  plane  strain  shear 
straining  of  a  composite  cylinder.  Determination  of  matrix  stresses  due 
to  applied  coay>o8lte  transverse  normal  stress  is  obtained  from  the  plane 
strain  isotropic  and  transverse  shear  straining  solutions  by  superposition 
in  the  following  manner  (see  Figure  39  and  reference  (g)):  First,  the  "upper 

_ I 

bound"  solution  for  a  transverse  shear  stress  of  022/2  (Figure  3®a)  is 
transformed  to  axes  at  45°  to  the  2-3  axes  giving  a  tensile  and  compressive 
plane  strain  biaxial  stress  state  (Figure  39b).  Next,  the  plane  isotropic 

_ r 

straining  solution  corresponding  to  a  stress  of  a 22 ^  (Figure  39c)  is 

added  to  the  transverse  shear  solution,  giving  the  plane  strain  transverse 

_ \ 

stress  solution  (Figure  39d) .  The  resultant  composite  axial  stress  0^ 

required  to  maintain  plane  strain  is  determined  by  integrating  0^  over 

the  X.  *  constant  end  surfaces  of  the  element,  and  the  solution  to  an 

*  ^  « 

axial  compressive  stress  of  0^  (Figure  39e)  is  superimposed  in  the  plane 
strain  transverse  stress  solution  to  yield  the  desired  plane  stress  trans¬ 
verse  stress  solution  (Figure  39f) .  Reference  (g)  gives  more  detail.  The 
final  solution  for  the  transverse  normal  stress  problem  is 

ar“  -  a22  [(D1+D3/r4-2D4/r2)cos  20  +  ^  +  Ij/r2]  (III-6a) 

°69  “  °22  C“(»1+2D2r2+D3/r4)cos  20  +  I2  -  Ij/r2]  (lll-6b) 
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K*  is  composite  layer  effective  bulk  modulus,  L  is  phase  elastic  Lam£ 
parameter,  and  the  are  constants  found  from  the  equations 


C^A^  -  dC  (I  II  -6m) 

where  and  d^  are  given  in  Table  4. 

(4)  Solution  for  General  In-Plane  Composite  Stress.  Equations 
(III-3),  (III -5),  (III -6),  and  Table  4  give  stresses  in  matrix  material 
of  a  unidirectional  composite  subjected  to  average  composite  axial  normal 
stress,  axial  shear  stress,  and  transverse  normal  stress,  respectively. 
Since  the  equations  are  linear,  they  may  be  superimposed  to  give  the  matrix 
stresses  due  to  any  general  combination  of  in-plane  conposite  stresses. 

The  final  result  will  be  of  the  form 
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where  Che  usual  Index  noCaCion  is  used;  i  and  j  can  Cake  on  values  of  r, 

8,  or  z;  a  and  8  are  1  or  2;  0  is  applied  composite  layer  stress;  om  is 
matrix  stress;  and  the  are  functions  of  position  coordinates  r  and  6, 

elastic  constants  of  fiber  and  matrix,  composite  layer  elastic  constants, 
and  phase  volume  fractions.  The  A^Qg  are  determined  readily  from 
equations  (III-3,  5,  and  6)  and  Table  4. 

3.  VISCOELASTIC  CYCLIC  HEAT  GENERATION  MODELS.  It  is  assumed  that 
all  mechanical  energy  dissipated  during  cyclic  stressing  of  a  visco¬ 
elastic  matrix  composite  material  goes  into  heat.  In  order  to  determine 
the  energy  dissipated  in  cyclic  loading,  and  hence  heat  generation, 
several  models  were  considered. 

a.  First,  a  purely  empirical  exponential  relation  for  energy 
dissipated  per  unit  cycle  of  stressing  was  considered: 

Dc  -  C oomwn  (III-7) 

where  Dc  is  energy  dissipated  per  unit  volume  of  composite,  oq  is  a 
measure  of  stress  amplitude,  u  is  cyclic  frequency,  and  c,  m,  n  are  em¬ 
pirical  constants  which  can  be  determined  from  experimental  data.  This 
specific  form  was  chosen  since  it  is  physically  obvious  that  Dc  must  be 
zero  whenever  aQ  or  u  go  to  zero;  hence,  m  and  n  must  both  be  positive. 

The  primary  drawback  of  this  model  is  that  it  does  not  predict  a  maximum 
level  of  energy  dissipated  per  unit  time  as  cyclic  frequency  goes  to 
Infinity.  Since  this  is  a  property  of  all  linear  viscoelastic  materials, 
the  exponential  form  (III-7)  was  judged  unfit  for  the  present  analysis. 

b.  An  analysis  of  a  general  linear  viscoelastic  material  under 
cyclic  stresses  was  then  made.  The  differential  viscoelastic  stress-strain 
relation  for  linear  viscoelastic  materials  is  of  the  form  (reference  (h)) 

Po  -  Qe  (III- 8a) 

where  a  and  e  are  stress  and  strain  measures,  and  P  and  Q  are  linear 
differential  operators  with  constant  coefficients 
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Here,  m  and  n  are  any  integer  numbers,  p^  and  q^  are  material  (visco¬ 
elastic)  constants,  and  t  is  time.  It  is  knovn  that  epoxy  matrix  materials 
behave  like  viscoelastic  solids  over  moderate  time  intervals.  The 
simplest  of  these  is  the  three-parameter  solid  modeled  by  a  Kelvin-Voigt 
parallel  spring-dashpot  element  in  series  with  another  spring,  for  which 
equation  (II 1-8)  becomes 

o  +  PjO  *  qQe  +  q^e  (III-9) 

where  the  overdot  denotes  time  differentiation.  Although  higher  order 
parameter  solids  were  considered,  the  three-parameter  solid  has  all  of  the 
basic  behavior  desired  in  a  viscoelastic  solid,  and,  because  of  its 
simplicity,  will  be  easier  to  work  with  empirically.  It  is  also  felt  that 
use  of  a  higher  order  model  would  result  in  "overkill"  and  could  not  be 
justified. 


c.  With  the  fundamental  form  (£11-9)  chosen,  it  is  necessary  to 
select  the  form  of  the  stress  and  strain  measures,  o  and  e.  Although  in 
general  most  polymeric  materials  exhibit  viscoelastic  volume  change,  this 
volumetric  viscoelastic  behavior  is  generally  considerably  less  than  that 
due  to  pure  distortion,  especially  for  low  stress  and  strain  levels. 
Accordingly,  the  stress  and  strain  measures  chosen  for  the  analysis  were 
the  stress  and  strain  deviators,  and  e^  ,  respectively,  given  by 
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where  and  are  strain  and  stress  tensors,  and  6^  is  the  Kronecker 
delta.  The  usual  subscript  notation  convention  is  employed.  The  assumed 
constitutive  relation  for  the  matrix  material  is,  therefore,  the  Isotropic 
form 
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In  order  to  obtain  energy  dissipated  per  unit  cycle,  it  is  assumed  that  the 

o 

applied  deviatoric  stresses  are  in  the  form  of  a  cyclic  component  S.,  super- 

* 

imposed  on  a  constant  level  S^: 

Sy  ■  S*j  +  S°j  sin  at  (III-12) 


-  ft  i  Vi  Arfitt*'1 
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The  resulting  strain  deviator  can  be  found  by  using  (III-12)  and  solving 
(III-ll)  to  give 

(III-13a) 


ey  »  AySin  wt  +  By  cos  +  eij 


where 


Ay  •  Sy  [1  +  ^  2^***^^^ 

By  ■  Sy  (p1q1-q2)tt>/(q12+q22w2) 


eij  "  Sij/ql 


(III-13b) 

(III-13c) 

(III-13d) 


Energy  dissipated  per  unit  cycle  per  unit  volume  of  matrix  D  is  given  by 

U>t»U>T  C 


Dc  “  JSijdeij  "  <*>  Sij4ij 


d(u>t) 


(III-14) 


cycle 


which,  upon  substitution  of  equations  (111-12,13)  into  (III-14),  becomes 


D«  _  n«(q2-piqi)  so  so 

c  “2 - 2~2  iJ  ij 

qx  +  q2  « 


(III-15) 


Equation  (III-15)  gives  the  energy  dissipated  per  unit  volume  in  one  cycle 
of  stressing  for  the  matrix  material.  The  composite  material,  treated  as 
an  "average"  material,  heats  up  by  an  amount  equal  to  the  energy  dissipated 
per  unit  composite  volume  multiplied  by  the  volume  of  composite  under  cyclic 
stress.  The  energy  dissipated  in  a  cycle  of  loading  per  unit  composite 
volume,  Dc,  can  be  found  by  averaging  the  matrix  energy  dissipated  over  the 
composite  volume  V: 


>c  - 1  i°c 


dV  - 


™(q2-P1ci1> 


2,  ,2 


qx  +  q2~“ 


TrcKq^Pjq^ 


2  . 

ql  +q2 


_  c°  s° 

V  ij  ij 


?Js«s« av 


(III-16) 


(III-16) 


where  the  overbar  denotes  volume  average.  The  energy  dissipated  per  unit 
composite  volume  per  unit  time  (energy  dissipation  rate  density)  is  found 
by  dividing  D  by  the  time  for  one  cycle,  2ir/u>,  to  give 


D  -  D  / (2ir/(ii)  - 


(q2*Plql>“ 

2  2~T 

2<ql  +  q2  “ 


S°  S° 

ij  ij 
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(III-17) 


d.  The  problem  of  determining  either  energy  dissipation  density  per 
cycle,  Dc,  or  energy  dissipation  rate  density,  D,  Is  seen  to  reduce  to  one 
of  determining  the  conposlte  volume  average  of  the  inner  product  of  the 
stress  deviator  amplitudes,  s°jS°j,  and  the  material  viscoelastic  param¬ 
eters  q^  and  p^.  As  mentioned  previously,  epoxy  materials  are  only 
weakly  viscoelastic  at  low  stresses  and  moderate  time  spans.  Therefore, 
an  extremely  good  approximation  to  the  viscoelastic  matrix  stress  ampli¬ 
tudes,  in  terms  of  applied  composite  layer  stress  amplitudes,  o°. 

ij  an 

will  be  given  by  the  elastic  stress  aeplitudes  found  from  equations 
(111-3,5,6,  and  7)  and  Table  A  of  the  preceding  section*  Equation  (111-10) 
is  then  used  to  determine  the  matrix  stress  deviator  amplitude  S°  . 


Finally,  the  composite  volume  averaged  is  determined.  Results  will 

be  of  the  form 


o  o 

SijSU 


o  o 

a&y6°a&ay& 


(III-18) 


where  the  constants  depend  upon  fiber  and  matrix  elastic  constants 
and  phase  geometry.  The  remaining  task,  that  of  finding  matrix  visco¬ 
elastic  properties,  is  handled  empirically  as  follows. 

e-  Equations  (III-15  and  17)  can  be  rewritten  in  the  form 


o  o 

it)  Q  C 

si-1si1 
A  +  Bw2 


“2 


2iKA+Bu  ) 


(III-16a) 


(III-17a) 


where 


A  -  qx  Mqj-Pjq^ 
2 

B  -  q2  /w(q2“P1<l1) 
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It  is  clear  that  If  Che  only  quantity  of  interest  is  heat  energy  generated 
and  not  the  general  viscoelastic  behavior,  then  only  the  two  constants  A 
and  B  need  be  determined  instead  of  the  three  p's  and  q's.  The  constants 
A  and  B  can  be  found  from  mechancial  cyclic  tests  on  specimens  of  a  given 
material  system  by  determining  the  energy  dissipated  per  cycle  from  the 
composite  stress-strain  curve  hysteresis  loop  at  several  cyclic  stress 
amplitudes  and  frequencies.  The  next  section  describes  such  cyclic  tests 
performed  on  graphlte/epory  and  S2-glass /epoxy . 

4.  CALIBRATION  TEST  PROGRAM  AND  MODEL  EVALUATION 

a.  The  empirical  viscoelastic  matrix  constants  A  and  B  in  the  energy 
dissipation  equations  (III-16a,  17a)  have  been  experimentally  determined 

for  the  graphite /epoxy  and  S2- glass/epoxy  composites  described  previously. 
Cyclic  tests  were  attempted  on  unidirectional  and  [±45]s  laminates  of  both 
materials.  Specimens  were  strain  gaged  and  tested  on  a  cyclic  tensile 
testing  apparatus  specially  built  for  this  purpose.  Load  cell  and  strain 
gage  outputs  were  amplified  and  fed  into  an  x-y  oscilloscope.  The  re¬ 
sulting  hysteresis  loops  were  photographed  and  energy  dissipation  density 
per  cycle  calculated  for  several  levels  of  applied  composite  load  and 
cyclic  frequency.  A  linear  regression  analysis  was  performed  on  the  data 
and  the  constants  A  and  B  determined. 

b.  Tensile  specimens  were  cut  from  panels  of  [0]  and  [±45]  laminates 
to  provide  axial  tensile  (fibers  in  load  direction),  transverse  tensile 
(fibers  perpendicular  to  load  direction)  and  primarily  axial  shear 

([±45]  angle  ply)  specimens.  It  was  found  early  on  that  the  axial  tensile 
specimens  did  not  provide  a  measurable  hysteresis  loop  and  were  dropped 
from  the  testing  program.  Difficulty  with  breaking  the  relatively  fragile 
transverse  specimens  eliminated  them  from  providing  meaningful  results. 

Nearly  all  results,  therefore,  were  obtained  from  8-ply  graphite/  and  4- 
ply  S2- glass/epoxy  [±45]s  tensile  specimens  in  which  the  layers  are  in  a 
stress  state  dominated  by  axial  shear.  All  specimens  were  aopropriately 
end-tabbed  with  glass/epoxy  strips.  Dimensions  of  the  specimens  are 
shown  in  Table  5. 

c.  The  cyclic  tensile  testing  apparatus  used  was  designed  and  built 
specifically  for  this  program.  The  apparatus  consists  of  a  350-kg  vi bra- 
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tion  isolation  platform  conatructad  of  bolted  plywood  and  railroad  ties 
and  supported  at  four  corners  by  helical  springs  (Figure  40).  The 
highest  natural  frequency  of  the  platform  and  machine  is  2  Hz.  A.  stiff 
support  structure  fabricated  from  welded  25-mm  steel  plate  is  bolted  to 
the  isolation  platform  and  holds  a  hardened  steel  fulcrum.  Resting  on 
the  fulcrum  is  a  hardened  steel  lever  arm,  one  end  of  which  is  attached 
to  the  specimen  grip  and  load  cell  hardware,  and  the  other  end  of  which 
is  pinned  to  a  load  actuating  rod.  The  load  actuating  rod  is  driven  by  a 
65-mm  dia  shaft  with  a  9-mm  eccentric.  The  drive  shaft  is  in  turn  driven 
through  a  belt-and  stepped- pulley  arrangement  by  either  a  1-hp  or 
a  3-hp  electric  motor.  (It  had  been  initially  planned  to  actuate  the 
lever  arm  by  an  electronic  shaker  with  a  75-kg  exciting  force  and  fre¬ 
quency  capability  in  excess  of  10,000  Hz,  but  untimely  failure  of  the 
shaker  forced  the  motor-and-pulley  design.)  The  lever  arm  position  is 
variable  and  provides  lever  arm  ratios  of  30:1  to  less  than  1:1  in  order 
to  control  the  applied  specimen  displacement  and  hence  load.  The  specimen 
gripping  arrangement  consists  of  orthogonal  pinned  joints  hung  from  the 
lever  arm  to  prevent  transmitting  bending  loads,  and  hardened  steel  bolted 
grips  at  the  specimen  upper  end.  The  specimen's  lower  end  is  bolted 
directly  to  two  strain-gaged  aluminum  plateB  which  are  in  turn  pinned  to 
a  threaded  collar  and  bolt  configuration  used  for  applying  preload  to  the 
specimen.  The  preload  mechanism  is  secured  to  a  32-mm  plate  which  is 
welded  to  the  bottom  of  the  fulcrum  support  structure  and  bolted  to  the 
isolation  platform. 

d.  Tests  on  graphlte/epoxy  and  S2-glass/epoxy  specimens  were 
performed  at  cyclic  frequencies  ranging  from  28  to  80  Hz  (an  underdesigned 
connecting  bolt  prevented  higher  frequencies)  and  specimen  average  composite 
tensile  stress  amplitudes  ranging  from  about  7.5  MPa  to  about  25  MPa. 
Prestressea  were  kept  at  about  1.5  times  the  amplitude  stress  to  avoid  com¬ 
pressive  loading.  Stresses  in  the  composite  were  kept  well  below  levels 
which  could  cause  matrix  cracking  in  the  specimens.  Each  specimen  was 
allowed  to  cycle  for  a  short  time  at  the  predetermined  stress  and  frequency 
to  allow  grip  setting  and  thermal  equilibrium  to  be  reached.  Load  cell 
output  and  specimen  strain  were  fed  to  the  ordinate  and  abscissa,  respec- 
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tively,  of  an  oscilloscope,  and  the  resulting  hysteresis  loop  photographed. 
Figure  41  shows  typical  hysteresis  loops  obtained.  Note  that  the  areas 
contained  in  the  loops  are  relatively  small.  The  area  contained  in  each 
hysteresis  loop  was  measured  using  a  planlmeter  and  the  mechanical  energy 
density  per  unit  cycle  was  calculated.  In  addition,  the  composite  volume- 
averaged  stress  deviator  amplitude  inner  product,  was  calculated 

from  equations  (111-3,5,6,7)  and  Table  4.  Fiber  and  matrix  phase  proper¬ 
ties  and  composite  properties  used  in  the  calculations  are  shown  in  Table 
5.  Test  and  calculation  results  are  shown  in  Table  6. 

e.  In  order  to  determine  material  viscoelastic  constants  A  and  B 
for  both  graphite/ epoxy  and  S2-glass/epoxy  material  systems,  the  energy 
dissipation  density  per  cycle  equation  (III-16a)  is  written  in  the  form 


(16b) 


The  quantity  S°^S°j/Dc  was  treated  as  the  dependent  variable  and  u>  as  the 
Independent  variable,  and  a  linear  regression  analysis  (reference  (i))  was 
performed  on  the  data.  The  result  is  a  determination  of  A  and  B  by  a 
least-squares  curve  fit  to  the  data  and  calculation  of  the  correlation 
coefficient,  an  estimate  of  the  "goodness  of  fit",  it  is  noted  that  the 
least-squares  curve  fit  is  unaffected  by  the  way  the  variables  are  struc¬ 
tured;  however,  strictly  speaking,  one  should  not  combine  dependent  with 
independent  variables  for  a  regression  analysis  as  was  done  with  D£  and 

Therefore,  the  correlation  coefficients  obtained  for  the  equation 
(III-16b)  are  not  valid,  and  are  not  presented  here. 

f.  Results  of  the  least-squares  curve  fit  to  the  data  and  resulting 
values  of  A  and  B  are  shown  in  Figure  42  for  graphite /epoxy  and  Figure  43 
for  S2-glass /epoxy .  Measured  values  of  0£  from  hysteresis  loops  are  plotted 
against  least-squares  curve  fit  predictions.  The  45°  line  represents 
perfect  correlation.  Note  that  there  is  considerable  scatter  to  the  data, 
most  likely  caused  by  the  difficulty  in  obtaining  accurate  planlmeter 
measurements  of  the  areas.  The  errors  are  estimated  to  be  as  high  as  +1002 
and  -502.  Even  with  this  high  probable  error,  the  correlation  is  good. 
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g.  With  tha  viscoelastic  constants  A  and  B  determined  from  the  ex¬ 
perimental  data,  an  analysis  can  be  performed  which  determines  the  maximum 
energy  dissipation  rate  density,  D,  for  a  given  value  of  applied  composite 
stress.  From  equation  (lll-17a),  it  can  be  seen  that  as  the  cyclic  fre¬ 
quency  a>  approaches  infinity,  the  energy  dissipation  rate  density  approaches 
an  upper  limit  given  by 

'  SyS^/2.11  (IH-19) 

Values  for  this  maximum  energy  dissipation  rate  density  are  shown  In  Table 

7  for  the  two  materials  tested.  Obviously,  there  will  be  an  optimum  finite 

frequency  where  most  of  D  will  be  obtained  (it  is  impossible  to  cycle  at 

max 

u  «  •) .  If  the  desired  dissipation  rate  density  is  arbitrarily  chosen  to 
* 

be  90X  of  the  cyclic  frequency  at  which  this  dissipation  rate 

density  will  occur  is  found  from 

•  «•»  -  “90  ^SV2*<A+B“«0> 

to  be 

u90  -  ICA/B)1* 

Values  of  0.9D  and  corresponding  cyclic  frequencies  are  shown  in  Table 
max 

7  compared  to  typical  values  obtained  from  a  test  machine  capable  of  35 -Hz 
frequency.  It  is  seen  that  the  higher  levels  of  energy  dissipation  rate 
density  are  much  larger  than  those  obtained  from  a  typical  cyclic  test 
machine.  This  has  implications  for  SGTF  thermographic  NDT  as  will  be  seen 
later. 

C.  EVALUATION  OF  CYCLICALLY  LOADED  GRAFHITE/EPOXY  AND  GLASS/EPOXY 

1.  SCOP!  OF  THE  EVALUATION.  It  was  necessary  to  determine  whether  or 
not  heat  generation  in  composites  due  to  viscoelastic  matrix  behavior  is 
adequate  for  thermal  imaging,  and  to  insure  that  the  heat  generated  near 
flaws  is  not  due  to  material  damage.  To  accomplish  this  task,  one  each  of 
2. 5-csMrt.de  x  150-cm  gage  length  samples  of  8-ply  [0/±45/90]  graphite /epoxy 
and  S2-glass/epoxy  were  drilled  with  a  6-mn  through-hole.  The  samples 
were  cycled  at  34  Hz  at  loads  below  those  which  would  cause  any  composite 
damage.  Thermo video graphs  were  taken  of  the  resulting  temperature  increase 
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of  Che  specimens.  Concurrently,  the  viscoelastic  heat  energy  generated 
was  calculated  for  both  specimens  using  the  equations  previously  presented 
and  the  general  procedure  outlined  in  Figure  36.  If  analysis  agrees  well 
with  experiment,  then  the  heat  generated  is  due  mainly  to  viscoelasticity 
and  not  to  material  damage.  If  experimental  specimen  surface  temperatures 
are  higher  than  predicted  near  the  hole,  then  material  damage  may  be 
creating  significant  heat,  thus  rendering  SGTF  methods  unsuitable  for 
Naval  aerostructure  thermographic  NDT. 

2.  FINITE  ELEMENT  STRESS  ANALYSIS.  The  8  ply  [0/±45/90]_  layup  is 
elastically  and  thermally  isotropic.  Elastic  and  thermal  properties  of 
the  [0/±45/90js  graphite/epoxy  and  S2-glass/epoxy  laminates  were  calculated 
using  laminate  analysis  techniques  (the  laminate  analysis  computer  program 
ELAN  was  used).  Table  8  presents  layer  and  laminate  thexmoelastic  proper¬ 
ties  for  the  materials  used.  Note  that  the  Poisson's  ratios  for  the  quasi¬ 
isotropic  graphite /epoxy  and  S2-glass/epoxy  are  the  same.  As  a  result, 

the  stress  distribution  in  the  test  specimens  containing  the  6-mn  through 
hole  will  be  similar  and  will  differ  only  by  a  constant  which  is  the  ratio 
of  the  applied  stress  on  the  two  specimens.  Only  one  specimen,  therefore, 
need  be  analyzed  for  stress  distribution.  The  tensile  specimen  containing 
the  hole  was  analyzed  by  the  SAFA  finite  element  computer  code.  A  relatively 
coarse  grid  was  utilized  since  the  object  is  to  obtain  temperature  dis¬ 
tributions,  not  peak  stresses.  The  finite  element  grid  representing  one 
quarter  of  the  specimen  due  to  double  symmetry  is  presented  in  Figure  44. 

The  analysis  was  performed  for  a  one  MPa  applied  laminate  far- field  stress, 
and  results  are  presented  in  Figure  45  averaged  over  adjacent  triangular 
elements.  Since  the  applied  laminate  stress  amplitudes  in  the  actual  tests 
were  13.4  MPa  for  graphite /epoxy  and  10.7  MPa  for  the  S2-glass/epoxy ,  the 
results  in  Figure  45  were  multiplied  by  13.4  and  10.7,  respectively,  to 
obtain  average  laminate  stresses. 

3.  LAMINATE  ANALYSIS  FOR  LAYER  STRESSES.  The  laminate  stress  analysis 
results  needed  to  be  analyzed  for  layer  stresses  in  order  that  the  heat 
generation/energy  dissipation  density  equations  could  be  used.  Therefore, 
a  constant  strain  laminate  analysis  was  performed  on  all  averaged  square 
element  stress  results  for  both  materials.  Again,  the  ELAN  computer  code 


:  -x. 
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was  used;  this  time,  to  obtain  composite  stresses  in  each  layer  of  the 
[0/±  45/90]  laminates. 

4.  HEAT  GENERATION  PREDICTION.  The  layer  stresses  were  used  in 
equations  111-3,5,6,7  and  Table  4  to  obtain  for  each  layer  of  each 

element.  The  computed  values  of  were  then  substituted,  with  values 

of  A  and  B  found  during  the  material  calibration  tests  and  a  cyclic  fre¬ 
quency  of  34  Hz  (corresponding  to  the  test  frequency) ,  into  the  energy 


dissipation  rate  density  equation  (III-17a).  The  results  were  heat  energy 
generation  rate  density  (equal  to  energy  dissipation  rate  density)  for 
each  layer  of  each  element  in  the  cyclically  loaded  laminate.  In  order  to 
provide  input  to  the  heat  transfer  analysis  used  to  predict  surface  tem¬ 
peratures,  heating  rate  density  was  averaged  over  layers  of  a  given  element. 
Results  are  shown  in  Figure  46  for  graphite /epoxy  and  S2-glass/epoxy . 

5.  SURFACE  TEMPERATURE  PREDICTION.  For  both  the  graphite /epoxy  and 
S2-glass /epoxy  samples,  a  51  node  model  was  created  for  one  quarter  of  the 
sample.  This  grid  system  was  possible  because  of  the  symmetry  in  the  sam¬ 
ple.  Figure  47  presents  the  nodal  pattern  and  significant  dimensions  of  the 
model.  Heat  generation  density  rates  from  Figure  46  were  used  with  the 
thermal  properties  in  Table  8  in  the  heat  transfer  analysis  to  predict 
surface  temperatures.  Results  show  that  although  the  S2-glass/epoxy  and 
graphlte/epoxy  samples  are  predicted  to  heat  0.9°c  and  0.1°C  above 
room  temperature,  respectively,  less  than  0.1°C  temperature  rise 
is  predicted  in  the  vicinity  of  the  hole  for  both  materials.  Thus,  a 
detectable  hot  spot  due  to  the  stress  concentration  effect  is  not  likely 
to  be  seen  in  the  thermovideographs  of  these  particular  cyclically  loaded 


specimens . 


6 .  EgPro^ENTAL^^RESlJhTS 

a.  The  graphlte/epoxy  and  S2-glass/epoxy  samples  containing  the 
6-tm  die  through-hole  were  cycled  for  30  minutes  each  at  34  Hz  and  tensile 
stress  amplitudes  ox  equal  to  13.1  MPa  andl0.7MPa,  respectively.  Figure 
48a  is  a  thermovideograph  of  the  S2-glass /epoxy  specimen  at  room  temperature. 
Black  dot  in  the  center  is  the  6r-mm  hole.  Note  that  the  specimen  tempera¬ 
ture  appears  about  0.1°C  cooler  around  the  hole  than  the  background. 

Figure  48b  is  a  thermo  videograph  of  the  same  specimen  after  cycling 
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for  30  minutes.  Two  aspects  of  the  figure  are  noteworthy.  First,  the 
lavender /red  region  in  the  center  is  about  0.6°C  higher  than  the  back¬ 
ground,  showing  that  the  specimen  as  a  whole  has  increased  about  0.7°C 
during  the  cyclic  loading.  Second,  there  is  no  measurable  temperature 
difference  between  the  area  immediately  around  the  hole  and  the  surrounding 
region.  The  graphite/ epoxy  tests  showed  similar  results  (Figures  48c,  d). 
Before  cycling,  the  region  surrounding  the  hole  appears  about  0.2°C  cooler 
than  the  control  strip  to  the  right;  after  cycling,  the  specimen  appears 
about  0.1°C  cooler  than  the  control.  This  indicates  that  the  specimen 
has  heated  about  0.1°C.  In  the  immediate  vicinity  of  the  hole,  there 
is  slight  evidence  of  a  hot  ring,  but  its  magnitude  is  less  than  0.1°C 
and  is  therefore  below  the  resolution  of  the  infrared  camera.  These 
results  agree  well  with  the  predictions  of  the  heat  transfer  analysis. 
Indicating  that  the  heat  being  generated  is  due  to  viscoelastic  behavior 
only,  and  not  to  continued  material  damage  due  to  cyclic  loading. 

b.  It  is  recognized  that  through-holes  and  partial  through-holes 
are  only  approximations  to  the  damage  which  typically  occurs  in  composite 
laminates  due  to  Impact  loading.  Actual  nonvisible  damage  (reference  (j)) 
consists  of  delamlnatlons  and  possible  fiber  breakage.  To  see  if  such 
damage  could  generate  hot  spots  visible  in  thermovideographs,  both  the 
S2-glass/epoxy  and  the  graphlte/epoxy  sample  were  impacted  by  a  blunt 
object  several  times  directly  around  the  hole.  Visible  delaminations 
resulted  in  the  hole  region.  The  samples  were  cycled  for  30  minutes  at 
the  previous  load  levels,  and  thermovideographs  taken.  Figures  49a  and  b 
present  the  results.  There  is  a  readily  visible  set  of  hot  spots  around 
the  hole  in  S2-glass/epoxy,  with  temperatures  as  high  as  0.2°C  above  the 
surrounding  area.  For  graphite/epoxy,  the  result  is  less  dramatic,  but 
present  nevertheless,  with  a  temperature  difference  of  about  0.1°C. 

7.  DISCUSSION 

a.  The  results  of  the  SGTF  study  show  that  at  cyclic  stress  levels 
below  those  which  can  cause  damage  in  composite  laminates,  the  heat 
generated  is  caused  by  viscoelastic  behavior  and  not  by  continuing  material 
damage  due  to  cyclic  loading.  However,  temperatures  were  not  significantly 
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higher  during  cyclic  loading  around  a  6-nm-dia  hole  drilled  in  25-mm- 
wide  tensile  specimens  of  S2-glass/ epoxy  and  graphit e/ epoxy .  This  is 
probably  due  to  the  relatively  low  levels  of  heat  generated  during  the 
tests,  and  would  most  likely  be  different  if  the  cyclic  frequency  were 
high  enough  to  produce  heat  energy  close  to  the  maximum  D  In  fact, 

as  seen  in  Table  7,  the  value  of  0. is  eleven  times  035  for  graphite/ 
epoxy,  and  more  than  three  times  larger  for  S2-glass/ epoxy.  It  remains 
to  be  seen  whether  this  heat  is  at  a  high  enough  level  to  produce  a  hot 
spot  where  there  is  no  delamination. 

b.  The  preceding  discussion  concerning  magnitude  of  heat 
generation  required  for  defects  not  containing  cracks  or  delaminations 
may  be  entirely  academic,  because  nearly  all  nonvisible  Impact  type 
damage  does  create  cracks  or  delaminations.  As  illustrated  by  the  SGTF 
tests  performed  on  impacted  specimens  with  holes,  the  surface  temperature 
near  holes  is  considerably  higher  than  without  delaminations.  The  reason 
for  the  higher  heat  generation  is  probably  frictional  heat  generated 
between  delaminated  and  cracked  surfaces.  This  is  reinforced  by  the 
preliminary  study  in  reference  (a)  which  showed  no  visible  heat  in  a 
slowly  cycled  unidirectional  specimen  of  graphlte/epoxy  containing  a 
6- mm  through-hole,  until  two  large  axial  cracks  developed.  There  was 
considerable  heat  generated  by  rubbing  of  the  cracks  long  after  the 
cracks  occurred.  (Figure  SO).  It  is  also  noted  that  Reif snider  and 
co-workers  (reference  (k))  have  been  able  to  generate  heat  in  delaminations 
by  applying  very  low  level  ultrasonic  excitation.  The  heat  occurs  only 
at  a  specific  (high)  frequency,  and  at  such  low  load  levels,  that  it 
cannot  be  viscoelastic  heat.  Friction  between  interfaces  is  therefore 
concluded  to  be  the  source  of  heat  for  their  tests. 

0.  FEASIBILITY  ASSESSMENT  OF  THE  SGTF  METHOD .  The  present  study  has 
shown  that  detectable  heat  is  generated  in  composite  materials  due  to 
viscoelastic  matrix  behavior  and  frictional  rubbing  of  interfaces  during 
cyclic  loading.  This  heat  is  generated  at  load  levels  which  should  not 
cause  material  damage  in  the  composites.  Therefore  the  SGTF  method  is  a 
safe  one  for  structural  composites.  Through-holes  without  delaminations 
cannot  be  detected  at  3 5- Hz  cyclic  frequency,  but  considerably  more  heat 
can  be  generated  at  a  few  hundred  Hz,  and  it  is  possible  that  such  defects 
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can  be  detected  at  the  higher  frequencies.  In  any  event,  composite  defects 
caused  by  impact  have  cracks  or  delaminations  associated  with  them;  and 
cracks  and  delaminations  can  produce  large  amounts  of  heat  due  to  relative 
motion  of  frictional  interfaces.  It  is  concluded  that  the  SGTF  method  is 
feasible,  but  determination  of  the  range  of  types  and  sizes  of  defects  is 
necessary  before  a  final  decision  can  be  made  concerning  its  advantages 
and  disadvantages  compared  to  the  EATF  method.  At  present,  the  EATF  appears 
to  have  the  advantage  of  simplicity  of  heat  source  and  application  tech¬ 
nique,  but  it  is  not  known  if  SGTF  methods  can  detect  flaws  which  are  not 
detectable  by  EATF  methods  and  vice-versa. 

E.  ULTRASONIC  STRESS  6ENERATI0N  METHODS 

1.  A  brief  investigation  was  conducted  into  the  use  of  ultrasonic 
exciters  as  a  stress  generated  heat  method.  There  are  two  ways  that  the 
ultrasonic  exciter  can  be  used:  as  a  high  frequency  stress  applicator, 
and  as  a  flaw  natural  frequency  exciter. 

2.  The  only  known  work  to  date  using  ultrasonic  devices  to  generate 
visible  heat  patterns  in  composites  has  been  performed  by  Henneke,  Reif- 
snide.-,  and  co-workers  at  VPI  (reference  (k)>.  They  report  the  use  of  an 
ultrasonic  exciter  to  which  is  attached  a  panel  of  composite  material  con¬ 
taining  delaminations.  When  certain  ultrasonic  frequencies  are  reached, 
heat  is  generated  at  the  delamination  site.  Thermovideographs  of  the  area 
clearly  show  the  presence  of  the  delamination  and  its  shape.  A  heat 
transfer  analysis  was  conducted  on  the  composite  assuming  that  the  delami¬ 
nation  Interface  is  a  heat  source.  Results  show  that  the  surface  tempera¬ 
tures  predicted  agree  well  with  the  thermovideographs.  This  outcome  is 
consistent  with  frictional  surface  heat  generation:  If  it  were  due  to 
viscoelastic  behavior,  only  the  edges  of  the  delamination  would  heat  up, 
not  the  whole  Interface.  In  addition,  the  heat  generated  would  increase 
to  a  maximum  at  several  hundred  to  1000  Hz,  then  remain  nearly  constant 
with  Increase  in  cyclic  frequency.  Since  the  heat  was  generated  at  only 
certain  high  frequencies,  it  is  probably  due  to  excitation  of  a  delamination 
natural  frequency. 

3.  Since  cyclic  stress  frequency  for  90Z  of  maximum  viscoelastic  heat 
generation  is  in  the  sonic  range,  going  to  ultrasonic  frequencies  for 
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viscoelastic  heat  generation  is  not  warranted.  Also,  considerably  larger 
energy  inputs  would  be  required  to  ultrasonlcally  generate  stresses  high 
enough  for  thermographically  visible  viscoelastic  heat  generation. 

4.  The  only  viable  way  to  use  ultrasonic  methods  for  SGTF  flaw 

i 

detection,  then,  is  to  excite  natural  frequencies  of  vibration  of  delami¬ 
nated  regions  in  composites.  Since  the  units  tested  at  VPI  were  used  on 
a  sample  not  greater  than  30  cm  square,  it  appears  that  a  unit  capable  of 
testing  large  areas  of  aerostructural  composite  surface  may  need  to  be  of 
considerably  greater  size.  Further  investigation  into  the  ultrasonic  SGTF 
technique  is  warranted  in  order  to  fully  determine  its  limitations. 
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IV.  SUMMARY  AND  CONCLUSIONS 

A.  The  present  Phase  I  study  was  critical,  as  it  decided  whether  or  not 
SGTF  techniques  can  be  used  for  thermographic  NDE,  made  a  definitive 
assessment  of  the  suitability  of  EATF  methods  to  detect  different  types 
of  flaws  in  the  most  important  aerostructural  materials,  and  evaluated 
the  effects  of  surface  reflectivity  and  emissivity. 

B.  EXTERNALLY  APPLIED  THERMAL  FIELD  (EATF)  DEVELOPMENT  AND  EVALUATION 

1.  The  developmental  effort  performed  in  Phase  I  was  to  determine  the 
size  and  location  of  flaws  that  can  be  detected  in  various  thicknesses  of 
structural  materials,  and  to  evaluate  different  methods  and  magnitudes  of 
external  heat  application.  In  addition,  limitations  on  infrared  camera 
sensitivities  were  examined;  and  the  effects  of  surface  reflectivity  and 
emissivity  caused  by  paints,  abrasions,  etc.,  were  evaluated.  The  work 
was  performed  using  a  closely  coordinated  experimental/analytical  approach 
where  a  test  program  was  conducted  on  selected  samples  of  flawed  and  un¬ 
flawed  graphite /epoxy,  glass/epoxy,  Kevlar  49/epoxy,  and  aluminum,  and 
correlated  with  a  heat  transfer/ thermal  analysis  that  can  be  run  with  any 
size  and  type  of  basic  structural  material  or  flaw. 

2.  The  test  program  consisted  of  applying  external  heat  to  test 
samples  containing  flam  not  visible  to  the  naked  eye,  and  observing  and 
photographically  recording  thermal  surface  temperature  perturbations  with 
the  Navy's  IGA  Thermovision  infrared  camera.  Material  thicknesses  were 
varied  up  to  5  mm,  and  flaw  types  Included  delaminations  and  partial  through- 
holes  which  simulate  manufacturing  and  impact- induced  defects.  Several 
types  of  heat  sources  were  evaluated.  Including  forced  hot  air,  radiation 
heating,  and  conductive  heating  with  thermal  blankets. 

3.  Conclusions  from  the  test  and  analysis  program  are  as  follows: 
a.  Simulated  delaminations  and  partial-through-holes  can  be  de¬ 
tected  In  sizes  as  small  as  6  mm  across  and  located  more  than  halfway  through 
32-ply  samples  of  graphlte/epoxy  and  16-ply  samples  of  glass/epoxy  mater¬ 
ials.  As  yet,  the  upper  limit  on  sample  thickness  has  not  been  found. 

There  appears  to  be  no  upper  limit  on  detectable  flaw  size,  but  there  is 
probably  a  lower  limit  on  flaw  size. 
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b.  Laminate  anisotropy  appears  to  be  unimportant  to  quality  and 
duration  of  thermal  Images:  virtually  identical  results  have  been 
obtained  with  radiative  heat  sources  on  [0/ ±4 5/90]  (thermally  isotropic) 
and  [Og/iAS]  (thermally  anisotropic)  graphite/epoxy  structural  composites. 

c.  EATF  heat  generation  using  variable  energy  radiant  heat  sources 
has  been  found  to  be  the  best  method  for  flaw  detection.  Use  of  thermal 
blankets  in  the  form  of  electrical  resistance  strip  heaters  produced 

poor  results.  It  has  been  concluded  that  the  so-called  "conduction" 
method  of  heat  application  (heating  a  remote  section  of  structure  and 
allowing  heat  to  conduct  along  the  surface  to  be  examined)  is  inappro¬ 
priate  for  detection  of  delaminations  and  partial-through-holes  (impact 
damage) . 

d.  Radiant  heat  source  application  times  of  up  to  one  minute, 
and  radiant  energy  of  up  to  one  watt  per  square  centimeter  may  be  re¬ 
quired  to  detect  flaws  in  typical  Navy  aerostructures  made  from  composites. 

e.  The  effects  of  Navy  aviation  structural  surfaces  on 
reliability  of  results  can  be  significant: 

(1)  Aluminum  surfaces  painted  with  polymeric  Navy  aviation 
paint  show  "hot  spots"  under  infrared  thermography  which  are  due  to 
uneven  paint  thickness,  and  which  could  be  misinterpreted  as  a  delamina¬ 
tion  or  flaw  in  a  graphite/aluminum  structure.  Painted  aluminum  matrix 
composite  structures,  therefore,  are  probably  not  good  candidates  for 
thermographic  NDE.  Paint  thickness  variations  do  not  cause  hot  spots 
when  applied  to  polymeric  matrix  composites  such  as  graphite/ epoxy  or 
Kevlar /epoxy. 

(2)  Tests  have  shown  that  reflection  of  a  spurious  heat 
source  from  a  naval  aviation  painted  surface  can  also  be  misinterpreted 
as  a  flaw  by  the  infrared  camera.  Structural  surfaces  to  be  examined 
must  therefore  be  shielded  from  sources  of  possible  heat  reflection  such 
as  sun,  light  bulbs,  and  engine  exhausts. 

(3)  Black,  orange,  and  white  Navy  aviation  paint  appear  to 
be  at  different  surface  temperatures,  probably  due  to  different  emissivi- 
ties.  Black  appears  significantly  hotter  than  orange  or  white.  Orange 
and  white  appear  very  close  in  temperature,  but  the  effect  is  still 
visible.  Paint  color  borders  may  therefore  be  visible  on  the  thermo¬ 
vision  screen  while  examining  surfaces  for  flaws.  However,  since  they  are 
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easily  seen  by  ch«  examiner,  they  ehould  provide  no  source  of  error  in 
the  examining  process.  Even  the  least  qualified  sailor  will  recognize 
that  flaws  do  not  come  in  the  shape  of  "MO  STEP". 

(4)  Abraded  or  soiled  surfaces  can  appear  as  cool  or  hot 
spots  even  in  the  absence  of  thermal  gradients.  These  effects  can  be 
eliminated  or  reduced  by  making  sure  that  surfaces  to  be  examined  are 
clean  and  unroughened  prior  to  examination. 

C.  STRESS  GENERATED  THERMAL  FIELD  (SGTF)  DEVELOPMENT  AND  EVALUATION 


1.  It  was  known  that  partial- through-holes,  notches,  and  delamina¬ 
tions  can  be  detected  by  SGTF  methods  in  glass/epoxy,  boron/ epoxy,  and 
graphite/ epoxy  structures.  However,  it  was  not  known  whether  composite 
structures  can  be  damaged  by  the  levels  of  cyclic  stress  application 
necessary  to  generate  heat  for  thermal  flaw  detection.  Further,  the  range 
of  stress  amplitude  versus  cyclic  frequency  necessary  to  get  a  visible 
stress-generated  thermal  image  for  different  types  of  materials  was  not 
known.  As  a  result.  Phase  I  developmental  activity  determined  the  levels 
of  stress  and  frequency  which  can  be  applied  to  structural  materials 
without  causing  structural  damage. 

2.  An  analysis  was  performed  which  determined  stresses  in  matrix 
material  for  a  given  set  of  in-plane  composite  laminate  stresses.  This 
analysis  was  used  with  a  semi-empirical  equation  derived  from  linear 
viscoelasticity  to  determine  energy  dissipated  in  the  composite  matrix 
material  for  different  values  of  applied  composite  stress  and  cyclic  fre¬ 
quency.  The  empirical  equation  was  calibrated  by  cyclic  frequency  tests 
run  on  unflawed  [±45]  specimens  of  graphite/epoxy  and  glass/epoxy.  Hys¬ 
teresis  loops  formed  during  the  cyclic  stress  tests  were  due  only  to 
matrix  viscoelastic  behavior,  and  not  to  damage  propagation  or  fracture. 

The  area  contained  in  these  hysteresis  loops  represents  the  energy  dis¬ 
sipated  as  heat  during  cycling,  and  was  used  to  determine  the  empirical 
constants  in  the  equation  of  heat  generation  as  a  function  of  applied  cyclic 
stress  level  and  cyclic  frequency. 

3.  With  the  experimentally  determined  equation  for  heat  generation  in 
composites,  it  was  possible  to  calculate  the  heat  generated  near  struc¬ 
tural  flaws  due  to  viscoelastic  behavior  only  (no  damage).  Cyclic  stress 
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tests  were  run  on  specimens  of  [Q/±45/9Q]  graphite/epoxy  and  glass/epoxy 
containing  6-mm  through-holes,  and  measured  surface  temperatures  were  com¬ 
pared  to  those  which  would  be  caused  by  viscoelastic  behavior  only  (no 
damage).  Cyclic  testa  were  also  run  on  graphite/epoxy  and  glass/epoxy 
samples  where  delaminations  had  been  Induced  around  the  holes  by 
Impacting  with  a  blunt  object. 

4.  Conclusions  drawn  from  the  analytical  and  experimental  studies 
are  as  follows: 

a.  In  all  tests  performed,  load  levels  were  maintained  less  than 
those  required  to  cause  breakage  of  fibers  or  cracking  of  matrix  material. 
The  measured  heat  agreed  reasonably  closely  with  the  calculated  values, 
indicating  that  viscoelastic  behavior  of  the  material,  and  not  material 
damage,  was  causing  the  heat  generated  in  the  composites. 

b.  Temperatures  around  the  holes,  however,  were  not  significantly 
higher  than  the  surrounding  regions  of  the  specimens.  This  is  probably 
due  to  the  relatively  low  levels  of  heat  generated  during  the  tests 
which  were  performed  at  a  cyclic  frequency  of  35  Hz.  Calculations  show 
that  three  to  ten  times  more  heat  would  be  generated  at  frequencies  of 
160  to  350  Hz  for  S2-glass/epoxy  and  graphite/epoxy,  respectively.  These 
higher  levels  of  heat  may  be  adequate  to  produce  a  thermographically 
visible  heat  pattern. 

c.  Considerably  higher  heat  is  generated  by  frictional  Interplay 
between  delamination  or  crack  surfaces  than  is  generated  by  material 
viscoelastic  response.  Since  nearly  all  impact-type  damage  results  in 
delaminations  and/or  matrix  cracks,  viscoelastic  heat  may  not  be  required 
to  generate  visible  heat  patterns. 

d.  SGTF  techniques  appear  viable  to  develop  detectable  "hot  spots" 
in  flawed  composites  without  degrading  material  or  reducing  structural 
lifetime;  however,  further  tests  on  actual  damage  sites  at  frequencies 

of  several  hundred  Hz  are  required  to  resolve  the  situation. 

e.  Investigations  into  the  use  of  ultrasonics  as  an  EATF  or  SGTF 
heat  generation  technique  indicate  that  ultrasound  methods  may  not  be 
fruitful  due  to  the  size  and  complexity  of  the  unit  necessary  to  generate 
sufficient  energy  to  cover  a  large  surface  area;  however,  further  study 
is  warranted. 
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C.  GENERAL 

1.  Du*  to  conductivity  of facts,  graphite/ epoxy  Is  a  more  difficult 
notorial  than  S2-glass/*poxy  in  which  to  dovolop  visible  thermal  patterns 
by  both  EATF  and  SGTF  techniques.  For  graphite/ epoxy  with  EATF  techniques, 
visibility  times  on  the  order  of  a  fraction  of  a  second  to  one  minute  have 
been  observed.  However,  flaws  can  be  detected  provided  that  the  operator 
utilizes  a  high  heat  Input  and  is  familiar  with  the  type  of  thermal  Images 
to  be  expected.  S2-glass/apoxy  has  longer  image  visibility  with  several 
seconds  to  more  than  a  minute  for  EATF  heating. 

2.  Infrared  camera  sensitivities  required  for  flaw  detection  by  EATF 
methods  range  from  2°C  full  scale  to  50°C  (detectable  temperature  rises 
of  0.2°C  to  S°C).  Sensitivities  for  SGTF  methods  range  from  1°C  to  5°C 
(detectable  temperature  rises  of  0.1°C  to  0.5°C). 
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V.  CURRENT  STATUS  AND  RECOMMENDATIONS 

It  has  been  demonstrated  that  EATF  radiation  heating  techniques 
can  detect  delaminations  and  partial- through-holes  in  structural  com¬ 
posites,  and  that  SGTF  heat  generation  may  produce  significant  thermal 
patterns  to  allow  flaw  detection  without  creating  material  damage  or 
reducing  coaponent  lifetime.  While  considerable  knowledge  has  been 
gained  concerning  capabilities,  limitations,  and  test  techniques  required 
for  EATF  thermographic  NDE,  it  still  remains  to  establish  limits  on  de¬ 
tection  capabilities  and  flaw  resolution  for  actual  defects  rather  than 
simulations;  and  to  investigate  the  effects  of  backup  structures  such  as 
honeycomb,  spars,  and/or  ribs.  For  the  SGTF  technique,  it  is  now  necessary 
to  determine  levels  of  stress  and  cyclic  frequency  required  to  generate 
detectable  thermal  patterns  near  impact  damage  sites  and  delaminations, 
and  to  investigate  methods  of  safely  applying  high  frequency  cyclic 
loading  to  aero structural  components.  Radiant  heat  source  and  infrared 
detector  equipment  guidelines  also  need  to  be  developed  further  to  enable 
equipment  development  if  thermographic  NDE  is  adopted  for  Navy  aviation 
ground  support  usage. 
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Table  1.  AGA  Thermovision  System  680/1026  Color  Screen 
Temperatures  and  Sensitivities 


Temperature  Increase  Above  Black  Setting,  deg.  C. 


Camera  Sensitivity 
(Full  Scale  Temp. 
Range) 

■ 

2 

5 

10 

20 

50 

100 

200 

500 

1000 

■■■■■■ 

■■■■■■ 

hh 

■■■■■MB 

■■■■■■ 

■■■■■■■ 

■■■■■■ 

mm 

Screen  Color: 

mm 

Black 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Dark  Blue 

0.1 

0.2 

0.5 

1 

2 

5 

20 

50 

100 

Blue  Green 

0.2 

0.4 

1.0 

2 

4 

10 

40 

100 

EE] 

Green 

0.3 

0.6 

1.5 

3 

6 

15 

60 

150 

H 

Light  Blue 

0.4 

0.8 

2.0 

4 

8 

20 

40 

80 

200 

400 

Purple 

0.5 

1.0 

2.5 

5 

10 

25 

50 

100 

250 

500 

Lavender 

0.6 

1.2 

3.0 

6 

12 

30 

60 

120 

300 

600 

Red 

0.7 

1.4 

3.5 

7 

14 

35 

70 

140 

350 

700 

Orange 

0.8 

1.6 

4.0 

8 

16 

40 

80 

160 

400 

800 

Yellow 

0.9 

1.8 

4.5 

9 

18 

45 

90 

180 

450 

900 

White 

1.0 

2.0 

5.0 

10 

20 

50 

100 

200 

500 

1000 

and 

and 

and 

and 

and 

and 

and 

and 

and 

and 

above 

above 

above 

above 

above 

above 

above 

above 

above 

above 
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Table  2. 

Aluminum  Surface  Finish  Samples 

PAINT 

THICKNESS 

SAMPLE 

COLOR 

(cm) 

A 

Black/White 

.0038 

B 

Gray/White 

.0038 

C 

Orange /White 

.0038 

D 

Unpainted/ Whit e 

.0051 

E 

White/White 

.0051 

F 

White/White 

.0051 

G 

White/ White 

.0064 

H 

White/White 

.0038 

I 

White/White 

.0038 

J 

White/White 

.0038 

K 

White/White 

.0038 

L 

White/White 

.0064 

58 
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Table  3.  Fiber  Composite  Surface  Finish  Samples 


PAINT 

SAMPLE  PAINT  COLORS  CONFIGURATION 


*Gr-A,  K-A  White/White 


Gr-B,  K-B  Black/ White 


Gr-C,  K-C  Orange/ White 


Gr-D,  K-D  White 


Gr-E,  K-E  White/Black 


Gr-F,  K-F  White/ Black 


*Cr~Unldlr actional  Graphite/Epoxy ,  K-Unidirectional  Kevlar/ Epoxy 
**2-Te©  Coats,  l*One  Coat,  W*White  Paint,  0>0range  Paint,  ►Black  Paint 
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Table  5.  Model  Calibration  Specimen  Dimensions  and  Material  Properties 


SPECIMEN  DIMENSIONS  (mm) 


FIBER  PROPERTIES: 

Axial  Young' s  Modulus  (GPa) 
Trans.  Young's  Modulus  (GPa) 
Axial  Shear  Modulus  (GPa) 
Axial  Poisson  Ratio  v^2 

Trans.  Poisson  Ratio  v>. 


Graphite /Epoxy 

S2-Gla88 /Epoxy 

200. 

200. 

150. 

150. 

19.1 

25.4 

1.32 

0.98 

234. 

83. 

13. 

83. 

25. 

22. 

0.4 

0.2 

0.3 

0.3 

3.7 

3.9 

0.35 

0.37 

MATRIX  PROPERTIES: 

Young' a  Modulus  (GPa) 
Poisson  Ratio 


COMPOSITE  PROPERTIES: 

[0]  Layer: 

Axial  Young's  Modulus  (GPa) 
Trans.  Young's  Modulus  (GPa) 
Axial  Shear  Modulus  (GPa) 
Axial  Poisson  Ratio 
Fiber  Volume  Fraction 


NAEC-92-157 


Table  6.  SGTF  Calibration  and  Calculation  Results  - 
[±45]s  Graphite/Epoxy  and  S2-Glass /Epoxy 


M 

Cyclic 

Tensile  Stress 

Energy  Dissipation  Rate/Cycle 

T 

Frequency 

Half  Amplitude 

S^Si1 

1—  1  1  1  1  1  1  1  1  1  ■ 

L 

oi,  rad  (Hz) 

V* HPa 

AJ  ij 

(MPa)2 

—  "3 - 

D£,  N-m/m  cycle 

wma^m 

176  (28) 

132 

123 

■on 

402 

428 

217  (34.5) 

mm 

216 

147 

. 

450 

628 

■nn 

734 

588 

g 

15.61 

■Em 

915 

1384 

eu 

w 

289  (46) 

5.43 

mm 

363 

184 

£ 

11.20 

1213 

892 

£ 

361  (57.5) 

5.6 

4  88 

243 

PfjTl  1 

30.0 

586 

1300 

6 

28.4 

1765 

1230 

115.0 

5006 

4990 

13.58 

35.0 

1839 

1775 

A  *  6.66  TPa/sec,  B  ■  12.8  MPa-sec 

176  (28) 

4.09 

3.9 

287 

316 

7.84 

14.5 

867 

1162 

217  (34.5) 

ra 

■H 

479 

1783 

£ 

FBI 

1842 

1184 

o 

12.28 

1  35.4 

3044 

3171 

M 

289  (46) 

4.78 

■sp 

723 

524 

CO 

9.55 

■TW 

3466 

2096 

a 

V 

361  (57.5) 

5.80 

7.9 

1023 

783 

N 

CO 

9.55 

21.4 

2372 

2126 

11.94 

33.5 

3705 

3322 

17.91 

75.3 

4281 

7473 

A  ■  1.74  TPa/Bec,  B  ■  14.61  MPa-sec 

*  Froa  Empirical  Equation  Dc  ■  S°  S°  u>/ (A+Bu>2) 


u[rad/sec],  S°jS°^[(Pa)2],  A[Pa/sec],  B[Pa-sec] 
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Table  7.  Heat  Energy  Dissipation  Properties  of  Graphite /Epoxy 
and  S2-Glass/Epoxy  Composites 


Material 

Max.  Energy 

Dissipation 

Rate  Density, 

• 

D 

max 

-8  2* 

10  “  S 

o 

N-m/m  sec 

°*9  Vx 
IQ'8  s2 

O 

N-m/m  sec 

Cyclic 

Frequency 

for 

0.9D 

max 

“90 

rad/sec (Hz) 

Energy 

Dissipation 

Rate  Density 
at  35  Hz 

*35 

ID'8  S2 
N-m/m^sec 

Graphite 

/Epoxy 

1.24 

1.12 

2162(344) 

0.11 

S2-Glass 

/Epoxy 

1.09 

0.98 

1034(165) 

0.31 

*S 


s°iis°iy  <Pa>‘ 
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Table  8.  Composite  Thermoelastic  Properties  Used  In  Calculating 
Surface  Temperatures  of  [0/+45/90]  Specimens  with  6  nm 
Through-Holes. 


Graphi te /Epoxy 

S2  Glass/Epoxy 

UNIDIRECTIONAL  COMPOSITE  PROPERTIES 

Elastic  Properties  -  See  Table  5 

Thermal  Properties: 

Axial  Thermal  Conduct.  (W/m°C) 

80. 

0.7 

Trans.  Thermal  Conduct.  (W/m°C) 

0.5 

0.3 

r0/±45/90]  LAMINATE  PROPERTIES 

Elastic  Properties: 

Young's  Modulus  (GPa) 

55. 

22. 

Poisson's  Ratio 

0.31 

0.31 

Thermal  Properties: 

In-plane  Thermal  Conduct.  (W/m°C) 

40. 

0.5 

Trans.  Thermal  Conduct.  (W/tn  C) 

0.5 

0.3 

Specific  Heat  (J/kg °C) 

920 

960 

Density  (kg/m^) 

1140 

1880 

SAMPLE  A 


Figure  1 


Thertncvideograph  of  Black/White  Navy  Aviation  Paint  on 
Aluminum  Plate.  Black  on  Left  and  White  On  Right  of  All 
Three  Photos,  (a)  and  (b)  Show  Emissivity  Difference?  of 
Black  and  White  Paint,  (c)  Shows  Light  Bulb  Reflection 
At  Upper  Right.  Note  Streaks  Due  to  Paint  Nonuniformity. 


Figure  3.  Thermovideographs  of  White  Faint  on  Aluminum. 

(a)  Effects  of  Faint  Thickness  Variations, 

(b)  Reflection  From  Light  Bulb. 


i 
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(B)  Kevlar  49/Epoxy  (Sample  K-A) 


Figure  4.  Thermovldeographe  of  Composite  Plates  with  Single  and  Double 
Coats  of  White  Navy  Aviation  Paint. 
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(a)  K-B,  Room  Temperature 
(Typical  of  all  R.T.  Samples) 


(b)  K-B,  Heated  60s., 
Cooled  30s. 


(c)  Gr-B,  Heated  30s., 
Cooled  308. 


Figure  5.  Thermovideographs  of  Single  Coats  of  Black  and  White  Navy 
Aviation  Paint  on  Composite  Plates. 


(a)  K-C,  Heated  40s- ,  Cooled  70s.  (b)  Gr-C,  Heated  30a. ,  Cooled  120s. 


(c)  K-F,  Heated  30a.,  Cooled  150s.  (d)  Gr-F,  Heated  30s.,  Cooled  25s. 


Figure  6.  Tharaovideograptas  of  Single  Coats  of  Orange  and  White  Navy 
Aviation  Faint  on  Composite  Plates. 
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Heated  30«.,  Cooled  90s.  Heated  35s.,  Cooled  240s. 

(b)  Gr-D,  Non-Uniformities  are  Due  to  Heat  Application 


Figure  7.  Thermovideographs  of  Single  Coat  of  White  Navy  Aviation 
Paint  on  Composite  Plates. 
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Figure 


(a)  Boom  Temperature 


(b)  Heated  60s.,  Cooled  30s.  (c)  Heated  60s.,  Cooled  420s. 


8.  Thermovideographs  of  Sample  Gr-D  With  Clean  Abrasions  on 
White  Paint. 
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(a)  Sample  Head-On  Showing  Anomaly  Locations 


(b)  Camera  at  15°  Angle  Shoving  Light  Reflection 
From  Glossy  Paint  Except  In  Abraded  Area  at 
Top  of  Sample 


Figure  10.  Gr-D  White  Painted  Sample  Containing  Clean  Abrasion, 
Abrasion  with  Heavy  Coat  of  Soiled  Grease,  and  Light 
Coat  of  Soiled  Grease  on  Glossy  Surface. 


amssmz 
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GRAPHITE  EPOXY  (.238  CM  THK) 
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IED  INPUTS 
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Figure  17.  EATF  Sample  Mounting  and  Test  Fixture. 
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w  lir  ftf.. 


8* 


SOURCE  TEMPERATURE  <DEG  C>  M0 


Figure  22.  Radiant  Heat  Source  Designed  and  Built  For  EA.TF  Studies. 


I  1/8  THROUGH 

31  cm 
l 


Figure  25.  Thermovideograph  of  8-Ply  Graphite/Epoxy  Panel  Containing 
Delaminations  1  Ply  and  A  Plies  Down  From  the  Surface. 


I  a; 


Thermovideographs  of 


5 


Panel 


Illlt 
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Figure  29.  Thermovideographs  of  Flaws  A2,  B2,  and  C2  Located  Half  Through 
32-Ply  Graphite /Epoxy  Panel,  (a)  and  (b)  Show  Flaws  C2  (Right) 
and  B2  (Left),  (c)  Flaws  B2  and  A2. 
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Figure  30.  Theraovideographs  of  Flaws  A1  and  Cl  Located  28  Flies  Down 
in  32-Ply  Graphite/Epoxy  Panel,  (a)  and  (b)  Are  Two 
Different  Attempts  At  Outlining  Both  Flaws. 


Figure  31.  Thenaovideographs  of  Flaws  B1  and  Cl  Located  28  Plies  Down 
From  Front  Surface  In  32-Ply  Graphite/Epoxy  Panel,  (a) 

10  deg.  C  Sensitivity,  (b)  5  deg.  C  Sensitivity. 
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Figure  32.  Deface  Mep  of  Anisotropic 
As  Seen  Froa  Blind  Side* 


(b)  11mm  65Z  (LEFT),  11mm  50Z  (RIGHT) 
Delamination  (CENTER) ,  2°C  Sens . 


(c)  Same  as  (b),  5°C  Sens 


Figure  33.  Thermovideographs  of  Blind  Side  of  Anisotropic  [02/±453  Gr/Ep 
Panel  With  Partial- Through-Holes  and  Delaminations.  Percents 
are  Z  Through  From  Back  Side. 
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1 

75 
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V 
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95 
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I 
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15Z  THRI 

( 

’“I 

1 

95 

' 

1 

MB 

J 

| - yu  aa  - q 

/W  —  -] 

■*"  11  110  aa 

Figure  34.  Defect  Map  of  16- Ply  [0/±45/90]  S2-Glaee /Epoxy  Panel 
As  Seen  Proa  Blind  Side . 
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(a)  20mm  152  (LEFT),  11mm  602  (RIGHT)  (b)  Same  as  (a)  but  COOLED  360s. 
25mm  DELAM,  1  Ply  Down  (LOVER  L.); 

HEATED  20s.,  COOLED  10s. 


(c)  25mm  DELAM,  HALF  THRU, 
HEATED  30s.,  COOLED  10s. 


(d)  LEFT-TO-RIGHT;  20mm  152, 
11mm  602,  11mm  302  PTH, 
HEATED  10s.,  COOLED  10s. 


Figure  35.  Thermovideo graph  of  Partial-Through- Holes  (PTH)  and  Delaminations 
(DELAM)  in  16-Ply  £0/±45/90]g  S2-Glass/Epoxy  Panel.  Percentages 
are  2  Through  From  Back  Surface. 
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Figure  36.  Procedure  For  Cel cult ting  The  Teaperature  Distribution  In  A 
Composite  Laainated  Structure  Under  Cyclic  Loading. 


(c)  COMPOSITE  CYLINDER  ASSEMBLAC2) 


(d)  COMPOSITE  CYLINDER  GEOMETRY 
AND  POLAR  COORDINATE  SYSTEM 


KEY  FOR  (a) ,  (b) ,  (c) 


FIBER 


MATRIX- 


-m 


Figure  37.  Deterministic  Phase  Geometries  For  Unidirectional  Fiber 
Composites. 


•  • 
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I 
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(a)  TRANSVERSE  SHEAR  (b)  XFORM  TO  AXES 

(UPPER  BOUND)  @  w/4 : 


(c)  ADD  PLANE  ISOTROPIC  STRAINING: 


(d)  GET  PLANE  TRANSV.  STRAINING: 


(c)  ADD  AXIAL  COMPR.  STRAINING: 


l 


(£)  TRANSVERSE  NORMAL  STRESS: 


i 

Figure  39.  Method  of  Determining  Phase  Stresses  Due  to  Layer 
(  Transverse  Normal  Stress  Using  Superposition  of 

Transverse  Shear,  Plane  Isotropic  Straining,  and 
Axial  Compressive  Straining. 


i 

1 
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28Hz,  ox° 


7 . 84MPa 


u  ■  57.5Hz,  crx° 


11. 94MPa 


(a)  [±45]  S2-Glass/Epoxy 


■  ■ 
J 


(it  -  28Hz,  a*0  -  9.50MPa  u>  -  57.5Hz,  ax°  -  12.22MPa 

(b)  [±45]  Graphite/Epoxy 

Figure  41.  Typical  Hysteresis  Loops  Obtained  During  Energy  Dissipation 
Calibration  Tests  of  Graphite/Epoxy  and  S 2-Glass /Epoxy 
Composites . 
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D 

c 

Measured 

N-m/n 


cycle 

1000 


100 


17.1 

17.1 

• 

17.1 

17.1 

15.2 

15.2 

17.1 

17.1 

17.1 

17.1 

» 

15.2 

20.7 

17.1 

17.1 

I 

13.8 

13.8 

20.7 

20.7 

18.2 

16.3 

14.1 

13.8 

20.7  20.9  18.5  13.8 


24.8  24.8  22.9  16.0  5.5 


24.0  27.8  27.0  17.9 


24.5  29.8  32.5 


23.4  32.5  60.1 


INSULATED 


INSULATED 

NODE  WITH  HEAT  SOURCE  AND  SURFACE  CONVECTION 
CONDUCTIVE  COUPLING 
CONVECTIVE  COUPLING 

Figure  47.  Nodal  Pattarn  Uaad  In  Numerical  Haat  Transfer  Analysis  of 
Cycled  Gr/Ep  and  S2-Cl/Ep  Samples. 


(a)  S2-Glass/Epoxy  -  Mo  Cycles 


(b)  S2-Glass /Epoxy  @  30  min. 


(a)  Graphite /Epoxy  -  Mo  Cycles 


(b)  Graphite/ Epoxy  @  30  min. 


Figure  48.  Thermovideographs  of  [Q/±45/90j  Composite  Specimens  With 
6- am  Through-Hole  Under  Cyclic  Tension. 


(a)  S2-Glaas/ Epoxy 


(b)  Graphita/Epoxy 


flgura  49.  Thatmovidaographa  of  £0/145/90]  Conpoalta  Specimana 

Containing  Through-Holaa  and  Dalaainations  Gndar  Cyclic 
Tana ion. 
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